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ABSTRACT 

We present ALMA 870 /rm (345 GHz) data for 49 high redshift (0.47< z <2.85), luminous (11.7 < 
log(Lboi/L 0 ) < 14.2) radio-powerful AGN, obtained to constrain cool dust emission from starbursts 
concurrent with highly obscured radiative-mode black hole (BH) accretion in massive galaxies which 
possess a small radio jet. The sample was selected from WISE with extremely steep (red) mid- 
infrared (MIR) colors and with compact radio emission from NVSS/FIRST. Twenty-six sources are 
detected at 870 /im, and we find that the sample has large mid- to far-infrared luminosity ratios 
consistent with a dominant and highly obscured quasar. The rest-frame 3 GHz radio powers are 24.7 < 
log(P 3 .oGHz/lU 77 z“®) < 27.3, and all sources are radio-intermediate or radio-loud. BH mass estimates 
are 7.7 < log(MBH/-^ 0 ) < 10-2. The rest frame 1-5 /rm SEDs are very similar to the “Hot DOGs” (Hot 
Dust Obscured Galaxies), and steeper (redder) than almost any other known extragalactic sources. 

ISM masses estimated for the ALMA detected sources are 9.9 < log(MisM/Al 0 ) <11.75 assuming 
a dust temperature of 30K. The cool dust emission is consistent with star formation rates (SFRs) 
reaching several thousand Mq yr“®, depending on the assumed dust temperature, however we cannot 
rule out the alternative that the AGN powers all the emission in some cases. Our best constrained 
source has radiative transfer solutions with ^ equal contributions from an obscured AGN and a young 
(10-15 Myr) compact starburst. 

Subject headings: galaxies: active - galaxies: evolution - galaxies: jets - radio continuum: galaxies - 
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1. INTRODUCTION 

Gentral questions concerning coeval galaxy and super- 
massive black hole (SMBH) evolution include the rel¬ 
ative timescales and mechanisms for stellar mass and 
BH mass building including the roles and duty cycles of 
(a) secular vs. merger-triggered mechanisms for driving 
material into the central regions; (b) “radiative-mode” 
vs. “jet-mode” BH accretion modes and rates; and (c) 
“quasar-mode” vs. “radio-mode” feedback mechanisms, 
all as a function of epoch, galaxy mass and galaxy en¬ 
vironment. Jet-mode AGN are low-excitation systems 
in which the AGN is powered by advection-dominated 
accretion flows onto the BH with a low accretion rate. 


as recently reviewed by Heckman & Best (2014). AGN 
in jet-mode are expected to have low radiative emission 
and to be energetically dominated by the jet kinetic out¬ 
flow. Jet-mode radio-AGN are thought to be highly ef¬ 
fective in maintaining galaxies free of new gas and star 
formation via “radio-mode” kinetic feedback: jet infla- 
tion of bubbles in the surrounding hot intergalactic gas 
( [Croto n et al.||2006[ |Cattaneo & Teyssier||2007| |Fabian 
I 2 U1L)|). Radiative-mode, or “quasar-mode ’, AGiN have 
nigher accretion rates from a thin accretion disk whose 
radiation powers the narrow and broad line regions, and 
which is fed through a surrounding dusty torus or “torus¬ 
like” structure. Quasar-mode accretion has a lower duty 
cycle than jet-mode, occurring when large amounts of 
material are available for accretion onto the supermas- 
sive BH. Quasar-mode AGN are capable of powering ef¬ 
ficient feedback into the host galaxy via thermal winds 
from the accretion disk, disrupting star formation and 
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ejecting gas. 

Powerful jets are also found in a ~10% of radiatively ef¬ 
ficient AGN, the Radio-Loud QSOs (RLQs) (also known 
as Broad Line Radio Galaxies, BLRGs, and High Exci¬ 
tation Radio Galaxies, HERGs), and therefore jets may 
also contribute to feedback activity in qu asar-mode AGN 
(Holt et al. 2008 Nesvadba et al. 12008 1. It is often as- 
sumed that the high radio power Eanaroff-Riley type H 
(FRH) jets ar e too collimated to impact the ISM within 
host galaxies (De Young 20101 . However, high resolu ¬ 
tion hydrodynamic mo dels by|\Vagner &: Bicknell (2011), 
Wagner et al.| (|2012|) and | Wagner et ai.| (|2013 have 
demonstrated that the impact ot powerful jets on the 
ISM within the central regions of AGN host galaxies is 
highly dependent on the density, the filling factor and 
the average size of the cool clouds within the ISM. High 
porosity leads to the inflation of a cocoon by the over¬ 
pressured jet, leadi ng to a quasi-spherical bu bble being 


driven into the IS M. Wagner & Bicknell (20111 and Wag¬ 
ner et al. (2012|) find that moderate to high power jets 
can accelerate dense ISM gas from a hundred to several 
thousand km s“^, with wide-angle flows, within 10 to 100 
Myr. 


Best & Heckman (2012) and Best et al. (2014) have 
shown that the two populations ot radio- AG jN, quasar- 
mode (or radiative-mode) and radio-mode (or jet-mode), 
are both found across all radio luminosities, and that the 
radiative-mode radio-AGN show much stronger evolution 
with cosmic time than jet-mode radio-AGN, with an or¬ 
der of magnitude increase in space density out to z~l. 
This evolution in space density is similar to that of the 
star formation rate density and the quasar luminosity 
function, consistent with the scenario that the radiative¬ 
mode AGN are controlled by episodic cold gas supply. 

Gas flows associated with gas-rich mergers are the 
likely source of episodic cold gas fueling, important for 
building galaxies and BHs by triggering both starbursts 
and AGN activity. Morphological signatures of mergers 
have been found in a large faction of powerful, z < 0.7, 
radio galaxies (RGs) which display spectroscopic sig¬ 


natur es of young stellar populations (Tadhunter et al. 
20111 ), and those with the youngest stellar populations 


(< oil Gyr) show a correlation with mid- to far-infrared 
(MFIR) and [O III] lumin osity, indicating the pre sence 
of a radiative-mode AGN. Tadhunter et al.| (2011) also 
find signihcant complexity in correlations between the 
triggering or re-triggering of jets, recent star formation 
and the merging of the dual nuclei, implying chaotic gas 
infall histories during merger events, while [Dicken et al.| 
(2011) find evidence for the strongest correlation between 
recent star formation and radio jets for the radio-AGN 
with the most compact jets. These results support the 
idea that small (young) radio jets play an important role 
in the evolution of massive galaxies and SMBHs during 
merger-driven high accretion rate phases. 

In this series of papers we address the impact of 
young, moderate to powerful, radio jets from luminous, 
radiatively-efficient, highly obscured, radio-AGN, on the 
disruption of the ISM and star formation in their hosts 
at redshifts near the peak of galaxy and BH building, 
z~l — 3. We also consider the possibility of ISM com¬ 
pression and starburst triggering by jet kinetic energy. 
By selecting systems with a high mid-infrared (MIR) lu¬ 
minosity we aim to identify radiatively efficient AGN, 


and by selecting compact radio sources we favor young 
radio jets which are confined within the hosts. By select¬ 
ing AGN which are very red through the optical-MIR 
we favor highly obscured systems likely to be in a peak 
fueling phase. 

A nearby example of such a system is the MIR-bright 
QSO Mrk 231, which has a luminous radiative-m ode 
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2003), and powerful 

molecular outi 
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Fischer et al.|20 
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nuclear geometry. 

1 .1. Evidence for Radio Jet Interactions with Molecular 

Gas 

Most early studies of outflows from AGN targeted the 
ionized gas, which, because of its much lower mass, re¬ 
quires much less energy to disperse than the neutral and 
molecular gas. More recent work shows that jet-induced 
feedback can indeed impact both warm and cold molec¬ 


ular gas in radio galaxies (Feruglio et al. 


2010 Fischer 


et al. |2010[ [Sturm et al.||2011[ [Uasyra Gombes| 


Gombes et al.||2013[ [Morganti et al. 12013 Garda-Burillo 


et al. 1120141 Dasyra et al. 

120141 IGonzalez- Alfonso et al. 

2014| I'i'adhunter et al. 

20 ' 

L4a|b 

Feruglio et al.||2015|), in- 

eluding turbulence and 

shock-excited H 2 emission in ra- 


dio galaxies (eg. Morganti et al. 2003; Nesvadba et al. 
2008,2011a,2011b; Guillard et al. 2012, 2015). Molec¬ 
ular hydrogen emission galaxies (MOHEGs) have large 
mid-IR H 2 luminosities which are too high for photo¬ 
dissociation regions, and which are mo st likely gener¬ 


ated by jet-generated shocks in the ISM ( Appleton et al. 
20061 [Ogle et ar]|2010[ |Guillard et al.|[20l2[ |Lanz et al. 


2015).^ Most of these radio galaxies are in a phase of 


radio-mode accretion without strong evidence for a con¬ 
current radiative mode AGN core. Quasar winds may 
also contribute to molecular outflows in radiative mode 
radio AGN (eg. Veilleux et al. 2013), such as those we 
study here, therefore our sample is ideal for studying the 
relative importance of these two feedback modes in early 
feedback phases of heavily obscured objects. 

1 .2. MIR-Selected Highly Obscured Quasars 

A rare class of highly obscured and luminous quasars at 
redshifts above 1 was revealed in follow-up studies of ex¬ 
tremely red (from t he optical to the MIR) sources found 
in Spitzer surveys (Lutz et al.|200'5 Yan et al.|20 05 
letta et al.|2006llDey et ai.|20U81|Lacy et al.|2ullD. 


Pol- 

_ hese 

systems can have luminosities over lO^'^L©, and deep X- 
ray observations have revealed Compton thick AGN in 
some (A^h > 1.5x10^"^ cm“^; Polletta et al. 2008). The 
reddest of these Spitzer-selecXed systems are so metimes 
ref erred to as Dust Obscured Galaxies, “DOGs” (Dey et| 


al. 


2008). 

"Tie Wide-field Infrared Survey Explorer (WISE), 
which covers the entire sky at 3.4, 4.6, 12, and 22 /rm 


bands) (Wright et al.||2010 

Jarrett et al. 

2011 

Cutri et 

al.|2012), has opened up the entire MIR s 
QSO searches by their MIR signatures dEi 

ky to obscured 

senhardt et al. 

20121 Stern et ai.|20I2l Wu et al.||2012 Asset et al.|2013t 

Yan et al.|20I3 Bridge et al.|2013 Jones et al.|20I4 Wu 

















































































































































































Jet Feedback in Obscured Radio-Loud Quasars 


3 


et al.|2014 Tsai et al.|2015 Stern et al.|2014 ). Although 
not as deep as the largest spitzer surveys {e.g. Lonsdale 
et al. 2004; Ashby 2009; Ashby 2013), WISE is sensitive 
enough to see the dust thermal emission from the most 
powerful quasars to redshifts >4. The first results from 
WISE follow-up of the reddest sources (selected without 
regard to radio brightness) have indeed revealed an ex¬ 
tremely IR-luminous population of high-redshift quasars, 
exceeding lO^^Lf?^ in bolometric luminosity in some cases 
(Eisenhardt et al. 2012 Wu et al. 2012 


ones et ar 11^01411 Wu et al.|:jO'i41|Assel ei al.M0i51 

Tsai et al.|2015). The ir bolometric lun iinosities are MIR- 


dominated, which led Wu et al. (2012) to dub t hem “Hot 
DO Gs” ( Hot Dust obscured Galaxies), while Bridge et 


al. (2013) investigate Lyo Blobs discovered around a high 


percentage of their red WISE sample, the WISE Lyman 
Alpha Blobs (WLABs). X-ray observations show only 
faint fluxes, consistent with highly obscur ed X-ray AGN 
(Stern et al.||2014 Pinconcelli et al.||2015). In this paper 
we will henceforth refer to the main dis covery papers for 
these WISE Hot DOGs and WLABs, ( Eisenhardt et al. 
2012 Wu et aL]|2012 Bridge et al.||2()13 ), as EWB12. 


1.3. This Work 

We present a snapshot survey of 49 sources with the 
Atacama Large Millimeter/submillimeter Array, ALMA, 
at 870 /im, which represents a southern sky subset of our 
sample of 156 radio-powerful (RP) obscured quasar can¬ 
didates. The sample has been selected to be unresolved 
in the NVSS and EIRS T radio surveys (jCondon et ah, 
1998| Becker et al.|1995 1 and ultra-red in the WISE MIr' 
survey, with similar selection criteria as the Hot DOGs 
(EWB12), to search for radio-jet dominated feedback 
from massive, obscured, quasars. We also present red- 
shifts obtained from optical and/or near-IR (NIR) spec¬ 
troscopy for 45 of these ALMA-observed quasars, and 
additional FlR-submillimeter photometry from other fa¬ 
cilities. 

The sample is described in Section and the obser¬ 
vations in Section The results are presented in Sec¬ 
tion followed by SED model fitting and derivation 
of luminosities and masses in Section [5l The discus¬ 
sion is in Section |6] and conclusions in Section [71 Seven 
sources with near-IR spectr oscopy from FIRE o n Magel- 
lan have been discussed by [Kim et al.j (|2013|). [Jones et 


al. (2015[) have published deep SCUBA 85U /rm imag¬ 


ing the the James Glerk Maxwell Telescope (JGMT) 
for 30 northern sources from our overall sample of 156 
RP quasars, detecting four, and finding an excess of 
serendipitously-d etected 850^^ sou rces in the fields on 
^1 Mpc scales. Silva et al. (2015) have found an ex¬ 
cess of serendipitously-detected 870/im sources in the 49 
ALMA imag es discussed here, in agreement with |Jones[ 
et al. (2015), although on smaller physical size scales of 
~150 kpc, perhaps indicating an excess of starbursting 
SMGs in the fields of these QSOs. Subsequent papers will 
address Jansky VLA 8-12 GHz imaging of the full sam¬ 
ple (Garol Lonsdale et ah, in preparation), and VLB A G- 
band imaging of 90 sources, including 33 from the ALMA 
sample (Colin Lonsdale et al., in preparation). NIR J 
and Kg imaging for a subset of the ALMA sample from 
VLT/ISAAC, and VLT/XShooter spectroscopy, will be 
presented in A. Blain et al. (in preparation). We adopt 
a cosmology with Hq = 71 km s“^ Mpc“^, Ha = 0.73 


and Hm = 0.27. 

2. SAMPLE SELECTION 

For our overall sample of ultra red, radio powerful, 
sources, we cross-matched the NVSS catalog, which cov¬ 
ers the sky north of J = —40°, with point sources from 
the WISE Allsky Catalog. The ALMA-observed sub¬ 
sample described here was selected at an earlier time 
when only the WISE Preliminary Catalog was available, 
as explained further below. Extended WISE sources 
were rejected based on the cataloged ext flag. We 
used positional information from the higher resolution 
but smaller area FIRST catalog where available, and 
excluded ±10° from the Galactic plane. Within this 
area of 28,443 square degrees of mutual NVSS-WISE 
coverage are 54,457 WISE sources which have a robust 
point source detection (SNR > 7) in the WISE 12 ^m 
and/or 22 /rm bands and an NVSS/FIRST 1.4 GHz 
match within a separation of < 7", the best compromise 
between completeness and reliability based on a random¬ 
ized match analysis. The entire sample is illustrated in 
the WISE 3.4 —4.6—12 /rm color-color diagram in Figure 
[2 where we plot the WISE-NVSS sample color-coded by 
922 = log(/ 22 /.tm// 20 cm), the 22 /xm q parameter, which is 
a measure of radio loudnes s (Appleton et al. 2004 Ibar et 
al.|2008 Ivison et al.|2007 ). The sequence of low-redshitt 


normal spirals and star bursts has blue {W1 — W2) colors 
(where WX is the Vega magnitude in WISE band number 
X) with a wide range of (W2 — W2>) colors, a nd a cloud 
of AGN i s seen with redder {W1 — W2) colors (Wright et 
2010 Jarrett et al.[[20lT| Yan et al.[[2013 ). 'I'he most 


al. 


radio loud systems are seen towards the left of the each 
of the normal galaxy sequence and the AGN cloud, con¬ 
sisting of the radio mode and the radiative mode AGN, 
respectively. For low redshift galaxies the point source 
fluxes used for this hgure may underestimate the total 
fluxes. 

Highly obscured luminous AGN are expected to be 
among the reddest sources in this figure. Previous au¬ 
thors have used color cuts in MIR color space and/or 
MIR-optical space to select the reddest extragalactic 
Spitzer and WISE sources (eg. Dey et al. 2006; 
EWB12). We have instead chosen to reduce poten¬ 
tial bias caused by color cuts, due to the complexity 
of the observed MIR spectral shape of these sources, 
which depends strongly on redshift due to the PAH emis¬ 
sion features and the 9.7 and 18 /rm silicate features. 
Therefore we have included all sources lying significantly 
redward of the main WISE populations. This is illus¬ 
trated in Figure by the dashed line, defined empiri¬ 
cally as: {W2 - W3) ± 1.25(1F1 - W2) > 7. The re¬ 
sulting number of sources in the sample redward of our 
MIR selection threshold in Figure [l] is 1858. To reject 
radio-quiet systems we added the requirement that the 
log(/22/im//20cm) < 0 as illustrated in Figurej2 For com¬ 
parison, the color selection criteria for the Hot DOGs of 
Eisenhardt et al. (2012) and Wu et al. (2012) are either 
(1F2 — 1F4)>8.2 or {W2 — VF3)>5 .3, termed by these 
authors the “WlW2drop” criteria. [Bridge et al. ([20131 
have used a slightly different color selection inethod: 
{W2 — VF3)>4.8. Both studies also use brightness or 
SNR threshold similar to ours, and neither has a radio 
flux density criterion. 

We inspected all candidates in WISE and DSS images. 
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and SDSS images where available, rejecting low redshift 
galaxies, artifacts and confused sources. To avoid reject¬ 
ing galaxies or quasars that could plausibly be within 
the redshift range of interest, 1 < z < 3, no specific 
optical magnitude or MIR/optical flux ratio cut was im¬ 
posed. Galaxies with a size or brightness large enough to 
place them at z < 0.5 were rejected, while high surface 
brightness compact sources were retained to brighter lim¬ 
its as potential high redshift quasars. We are obtaining 
spectra to reject low redshift sources from our final sam¬ 
ple. Thus, our selection can include objects that are less 
red (in R-[22]) than the criterion {R — [24]) > 14 used 


for Spitzer DOGs (Dey et al. |2008) and similar sam¬ 
ples, and, for example, may include systems in which a 
massive galaxy dominates the optical light. A total of 
708 of the NVSS-WISE sample satisfy all our criteria, or 
^1.3% of the entire NVSS-WISE matched sample. Of 
the 708 sources, 703 (269) are detected at 12 (22) ^m 
with SNR>7. 



W2-W3 



Figure 1. WISE 3.4 — 4.6 — 12 (ira (W1 — W2 — W3) color-color 
plot (using Vega magnitudes) showing the full WISE-NVSS sample 
of 54,457 sources, color-coded by radio-loudness (see Figure]^ as 
shown by the color bar on the right, radio loudness increasing red 
to black. The sequence at the bottom is the sequence of normal 
spirals and star bursts. The flux-limited ultra-red sample of 156 
sources is highlighted above the dashed line with larger symbols. 
The hor izontal dashed line shows the AGN color-selection criterion 
used by|Stern et al.| (|2012||: (W1 — W2) > 0.8, for comparison to 
our ultra-red selection criterion. 


Due to the strongly varying WISE sensitivity across 
the sky, caused by the varying coverage level, a flux den¬ 
sity threshold at 22/tm was applied. Eor the subsample 
observed with ALMA the flux density threshold chosen 
was 4 mjy. This was later revised upward to 7mJy for 
the remainder of the sample when the Allsky Catalog 
became available. The final sample is 156 sources, 49 of 
which belong to the ALMA-observed subsample. 

The ALMA subsample fluxes have been revised in the 
updated WISE Allsky Catalog, some now falling below 
the original 4 mJy threshold. These have been retained 
in the sample. The declination range for the ALMA sub¬ 
sample was chosen to allow access to many northern facil¬ 
ities as well as ALMA: —40° < 6 <4-1°. The ALMA sam- 



log([22|j,m, mJy]) 

Figure 2. The 22/rm, q value, 522 = log(/ 22 Aim// 20 cm) vs. 
WISE 22 fim flux density, for the 156 ultra-red, radio-powerful 
sources (blue points), compared to the full WISE/NVSS cross- 
matched sample of 54,457 sources (dots and contours). The com¬ 
pleteness limit of the NVSS is responsible for the diagonal cut¬ 
off at the lower left. Also shown are the mean expected val¬ 
ues for star forming galaxies of the infrared-to-radio par ameter. 
Q9A = log ( f9Ai,.m / fcjOrin) for the SpUzev 24 p,m band by |Ibar et| 
\&i] 1 |2008[ | in the local universe and the k-corrected ^'^ 2. 0 value. 
Tne representative radio-loud SED from|Elvis et al.|l|1994||, is also 
shown. The subset of ultra-red sources observed witn ALiMA was 
selected to have 0 < log(/22^m//22cm) < -1- 


pie is restricted in RA to two regions: 3^ < RA < 8^30™ 
and 13** < RA < 21^, due to the incompleteness of 
the WISE Preliminary Catalog at that time. Fifty-five 
sources met our original selection criteria within these 
areas. Six of these have log(/ 22 pm// 20 cm) <-l and were 
excluded to disfavor classical double lobed sources, a cri¬ 
terion which was later dropped for the full sample of 156 
sources. Of the final ALMA subsample of 49 sources, 48 
(23) have SNR>7 at 12 (22) ^m in the WISE Allsky Cat¬ 
alog, and the minimum 12/22/rm SNR is 6.7/2. 6 . Only 
one source fails the original SNR criterion of SNR>7 
at 12 or 22 pm after the Allsky flux revision: WISE 
J204049.51-390400.5 which has an SNR at 12/22/im of 
Q.71&.2. Due to the evolving selection criteria the sample 
is not complete to a fixed SNR or flux density limit. 

3. OBSERVATIONS 

3.1. FIR and Submillimeter Observations 

Twenty-three sources were observed with ALMA in 
two Band 7 scheduling blocks on 2011 November 16. An 
additional 14 were observed on 2012 May 25 and the fi¬ 
nal 12 on 2012 August 29. The central frequency was set 
to 345 GHz (870 pm), with an 8 GHz bandwidth split 
into two sidebands. Fifteen antennas were used in the 
“compact” Early Science configuration for the November 
2011 observations, 19 for the May 2012 observations and 
23 for the August observations. The resulting beamsizes 
are 172, O'.'5 and 1'.'2, respectively. 

For the November observations, Callisto was used for 
flux calibration. Titan was used for 7 sources in May and 
Neptune for the other 7 sources in May. Titan was again 
used for all sources in August. The objects were ob- 
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served with the correlator in Time Division Mode, which 
results in 256 channels per sideband with a spectral res¬ 
olution of 14 km s“^. Time on source was 1.5 min per 
object, resulting in an rms noise of 0.3-0.6 mJy, the lower 
noise levels generally correspond to the observations with 
larger antenna numbers. Flux densities were measured 
using imfit in CASA, and are reported in Table 

Four sources from our overall sample of 156 were ob¬ 
served at 350 /im at the Submillimeter High Angular 
Resolution Camera II (SHARC-II) installed at the 10.4 
m Caltech Submill imeter Observatory (CSO) telescope 
(Dowell et al.|2003 ) on UT 2012 March 22. Two of these, 
W1343-1136 and W1400-2919, belong to the ALMA sub¬ 
sample. Our observin g and data reduc tion process fol¬ 
lows that described byfWu et al.|(|2012|) . 

Seventy-seven of our full sample of 156 sources, includ¬ 
ing all the ALMA targets, were queued for Hersche^\ 
PACS (Poglitsch et al. 2010) and SPIRE (Griffin et ^ 
2010) photometric observations (PID: OT2_clonsdaLl). 
In total 15 objects were observed by PACS, and 5 with 
SPIRE, before the cryogen depletion of the Herschel 
Telescope in 2013 April. For the PACS observations, two 
concatenated mini-scans orientated at 70° and 110° were 
acquired. Each mini-scan has 8 legs, a scan length of 3' 
and and a leg separation of 5", with a time on source of 
220 seconds. The SPIRE observations were conducted 
using the small jiggle map mode with 2 repeats, with 
74 seconds on source time. The data were processed 
with HIPE 11.1.0 (Ott 2010). The aperture-corrected 
flux densities for all 15 Herschel-observed sources are pre¬ 
sented in Table [21 

Summarizing, the FIR/submillimeter data available 
for the 49 ALMA-observed sources discussed in this pa¬ 
per, PACS observations are available for two sources, 
W1500-0649 and W2059-3541, and SPIRE observations 
for one, W2059-3541. Two additional sources have upper 
limits from the CSO at 350/rm, W1343-1136 and W1400- 
2919. 


3.2. Optical and Near-IR Photometry and Spectroscopy 

Redshifts for 45 of the 49 ALMA sources and i?-band 
photometry for 26 were obtained and are presented in Ta¬ 
bles |T] and 121 In total 48 of the 49 sources were observed 
spectroscopically, and three of these yielded no line detec¬ 
tions. Full details of the spectroscopy will be presented 
in a later publication; here we make use only of the red- 
shifts. Optical spectra of 23 objects were obtained using 
the Goodman spectrograph on the SOAR 4.2 m telescope 
on UT 2012 January 21-24 and UT 2012 December 9-12. 
The R-band photometry obtained with SOAR was taken 
with the spectrograph acquisition camera. We observed 
five sources from the WISE-NVSS-ALMA sample, and 
two from the northern JGMT sample, using the Double 
Spectrograph on the 5-m Hale telescope at Palomar Ob¬ 
servatory between 2012 Nov ember and 2013 M arch. As 
described in more detail by Kim et al. (20131, twenty- 
four sources were observed m the INIR with Magellan 
on UT 2012 July 27-29. Finally, twenty-eight sources 
were observed with VLT / XShooter on UT 2013 June 
04-06 and thirty-one sources were observed in Kg with 


Herschel is an ESA space observatory with science instru¬ 
ments provided by European-led Principal Investigator consortia 
and with important participation from NASA. 


VLT/ISAAC over 3 nights from UT 2013 June 01-04, of 
which 14 were also observed in J. 

4. OBSERVATIONAL RESULTS 

4.1. Photometry and Redshifts 

The R h r magnitudes and the WISE, NVSS and 
ALMA flux densities for the 49 sources are presented in 
Table HI There are 26 ALMA detections at 3cr or above. 
None of the sources are resolved. The two sources from 
the ALMA subsample observed with Herschel were both 
detected at 70/tm, and W1500-0649 also at 170/tm (Table 
Twelve additional sources from the full sample have 
Herschel detections. Neither of the two ALMA sources 
observed with the CSO, W1343-1136 and W1400-2919, 
were detected, resulting in 3(T upper limits of 45 mJy for 
each, while one of the other two souces, W1025-I-6128, 
has a modest detection, as listed in Table [2] There are 
available R-band Vega system magnitudes from SOAR 
for 16 of the ALMA-detected sources and for 10 of the 
sources with upper limits. W0823-0624 and W1434-0235 
have Sloan Digital Sky Survey (SDSS) data in one or 
more bands. The R/r magnitudes range from 20.7 to 
24.1. The optical photometry is used in this paper only 
as a rough constraint on the mass of the stellar popula¬ 
tions; a full analysis of these data and the NIR photome¬ 
try will be presented by A. Blain et al. (in preparation). 

Redshifts are available for 25 of the 26 detected 
sources, and for 20 of the 23 sources with ALMA up¬ 
per limits. Six of the 19 sources in Table have known 
redshifts, including the four sources in the ALMA sub¬ 
sample and the two JCMT-observed sources. It is be¬ 
yond the scope of the present paper to analyze the opti- 
cal/NIR spectroscopy, however we note that many spec¬ 
tra have indications from the ionization levels of an ob¬ 


scured radiative-mode AGN (Baldwin, Phillips & Ter- 
TMT] |Kewley|[2M3|. 

3ns have been publii 


I'he i VIagellan/llRE ob- 
published by Kim et al. (20131 


levich 
servations 

while the VLT/Xshooter observations will be reported 
by A. Blain et al. (in preparation). One object 
was not observed (W1657-0102), one was not detected 
even in continuum (W0614-0936), and for two objects 
only a faint continuum was detected (W1707-0811 and 
W2040-3541). Additionally the redshifts are uncertain 
for five sources which have only weak lines, a single 
line detection or an unresolved line pair. They are: 
W0519-0813 (LaA1216A); W0702-2808 ([O II]A3727A); 
W1521-L0017(weak H^,[0 III]A4959,5007A); W1703- 
0517 (blended Ha-L[N II]A6584A); and W1936-3354 
([O II]A3727A). 


4.2. Colors and Spectral Energy Distributions 

The ALMA observations were designed to constrain 
the luminosity of these quasars in the rest-frame FIR- 
submillimeter compared to the MIR. We will show that 
most of our WISE-NVSS-ALMA sample have SEDs dom¬ 
inated by an AGN in the MIR, and possibly through the 
FIR-submm also. However, substantial rates of star for¬ 
mation are likely also present. 

We present the main results in Figures [3j]^ We show 
these results in two complementary ways: (1) the ra¬ 
tio of 870 fim/22 /rm flux density as a function of red- 
shift is shown in Figure and (2) we use the range in 
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Table 1 

Optical, WISE, NVSS and ALMA Photometry 


WISE Name 

i2-band 
app. mag. 
(Vega) 

/3.4/im 

(mjy) 

WISE 

/4.6/jm 

(mJy) 

WISE 

fl2fj,m 

(mJy) 

WISE 

f22fiTa 

(mJy) 

WISE 

(122 

observed 

frame 

/l.4GHz 

(mJy) 

NVSS 

fsjOfim 

(mJy) 

ALMA 

log 

,/870/im / 
f22fi,ra) 

J030427.53-310838 

22.0±0.3 

0.179±0.0066 

0.693±0.018 

4.55±0.12 

10.87±0.65 

-0.70 

54.99±1.72 

2.8±0.3 

-0.52 

J030629.21-335332 

20.7±0.2 

0.194±0.0070 

0.749±0.020 

6.19±0.14 

4.83±0.65 

-0.11 

6.17±0.52 

1.7±0.4 

-0.45 

J035448.24-330827 

22.5±0.2 

0.086±0.0063 

0.314±0.015 

3.34±0.12 

7.05±0.78 

-0.04 

7.73±0.54 

<1.2 

<-1.12 

J040403.61-243600 

20.7±0.2 

0.104±0.0072 

0.326±0.017 

4.63±0.15 

4.43±0.92 

-0.40 

11.04±0.58 

3.1±0.4 

-0.16 

J040937.67-183757 

21.3±0.2 

0.517±0.0147 

2.492±0.054 

22.97±0.36 

47.07±1.68 

0.00 

46.91±1.48 

<1.5 

<-1.47 

J041754.10-281654 

23.9±0.3 

0.043±0.0057 

0.143±0.013 

1.43±0.10 

3.84±0.85 

-0.63 

16.32±0.68 

4.3±0.6 

0.05 

J043921.92-315908 

24.2±0.3 

0.034±0.0051 

0.090±0.011 

1.19±0.09 

4.57±0.77 

-0.66 

20.88±0.78 

6 ±0.4 

0.12 

J051905.84-081320 

>24 

0.059±0.0071 

0.136±0.015 

1.76±0.13 

5.09±0.85 

-0.72 

26.79±0.90 

<1.5 

<-0.56 

J052533.47-361440 

22.9±0.3 

0.013±0.0053 

0.046±0.010 

1.38±0.10 

4.07±0.71 

-0.10 

5.12±0.56 

<1.5 

<-0.41 

J052624.72-322500 

22.6±0.2 

0.023±0.0056 

0.106±0.011 

6.73±0.17 

26.47±1.17 

-0.82 

173.5±5.23 

18 ±0.5 

-0.17 

J053622.59-270300 

23.2±0.3 

0.065±0.0060 

0.404±0.016 

3.26±0.12 

5.35±0.80 

-0.18 

8.16±0.56 

2.7±0.4 

-0.30 

J054930.07-373939 

>24 

0.021±0.0049 

0.091±0.010 

1.09±0.09 

3.37±0.73 

-0.50 

10.60±0.57 

2 ±0.4 

-0.23 

J061200.23-062209 

20.8±0.2 

0.330±0.0118 

0.621±0.022 

9.52±0.20 

20.55±1.01 

-0.20 

32.79±1.08 

2.7±0.6 

-0.88 

J061348.08-340728 

>24.6 

0.065±0.0065 

0.211±0.012 

2.62±0.11 

7.37±0.71 

-0.52 

24.15±1.16 

<1.8 

<-0.62 

J061405.57-093658 

24.1±0.3 

0.017±0.0073 

0.069±0.015 

2.16±0.14 

5.80±0.84 

-0.01 

5.94±0.51 

<1.8 

<-0.51 

J063027.81-212058 

22.2±0.2 

0.020±0.0070 

0.063±0.014 

2.62±0.13 

5.06±0.90 

-0.39 

12.55±0.60 

5 ±1.3 

-0.01 

J064228.92-272801 

21.4±0.2 

0.038±0.0070 

0.108±0.014 

1.33±0.16 

3.67±0.87 

-0.24 

6.36±0.52 

2.2±0.6 

-0.22 

J065215.85-200612 

>23.7 

<0.013 

0.045±0.013 

1.92±0.13 

4.81±0.83 

-0.18 

7.35±0.52 

3.2±0.5 

0.18 

J070257.20-280842 

21.8±0.2 

0.027±0.0061 

0.140±0.016 

1.76±0.12 

4.54±0.86 

-0.25 

8.00±0.58 

<1.8 

<-0.40 

J071433.54-363552 

22.8±0.2 

<0.012 

0.039±0.011 

0.99±0.12 

4.01±0.84 

-0.47 

11.95±0.60 

2.4±0.3 

-0.22 

J071912.58-334944 

24.1±0.3 

<0.011 

0.081±0.012 

1.93±0.12 

4.06±0.88 

-0.78 

24.30±0.87 

5.2±0.6 

0.11 

J081131.61-222522 

21.5±0.2 

0.132±0.0086 

0.611±0.023 

5.62±0.17 

7.61±1.17 

-0.37 

17.84±0.71 

<1.8 

<-0.63 

J082311.24-062408 

22.32±0.16‘" 

0.118±0.0078 

0.441±0.019 

4.08±0.15 

10.42±0.97 

-0.71 

53.99±1.67 

<1.8 

<-0.76 

J130817.00-344754 

22.3±0.2 

0.086±0.0056 

0.248±0.013 

3.36±0.12 

9.12±0.73 

-0.87 

68.10±2.10 

1.38±0.34 

-0.81 

J134331.37-113609 

21.7±0.2 

0.024±0.0057 

0.136±0.013 

1.61±0.12 

3.81±0.79 

-0.33 

8.18±0.54 

2.34±0.31 

-0.22 

J140050.13-291924 

21.7±0.2 

0.110±0.0063 

0.501±0.018 

5.58±0.14 

11.85±0.77 

-0.64 

51.92±1.63 

<0.90 

<-1.37 

J141243.15-202011 


0.092±0.0063 

0.333±0.015 

3.39±0.13 

7.41±0.78 

-0.09 

9.01±0.55 

2.55±0.63 

-0.45 

J143419.59-023543 

22.02±0.18‘‘ 

0.058±0.0056 

0.257±0.014 

2.13±0.11 

5.04±0.71 

-0.86 

36.15±1.16 

<0.9 

<-0.75 

J143931.76-372523 


0.027±0.0071 

0.115±0.013 

2.34±0.12 

3.92±0.83 

-0.41 

10.02±0.57 

<0.6 

<-0.75 

J150048.73-064939 


0.068±0.0065 

0.293±0.016 

6.26±0.17 

15.77±0.94 

-0.10 

20.01±0.73 

6.11±0.28 

-0.41 

J151003.71-220311 


0.143±0.0095 

0.411±0.020 

5.34±0.18 

14.87±1.09 

-0.06 

17.14±0.70 

<0.9 

<-1.27 

J151310.42-221004 


0.037±0.0082 

0.214±0.018 

2.64±0.16 

9.71±1.10 

-0.50 

30.40±1.03 

4.86±0.27 

-0.30 

J151424.12-341100 


0.076±0.0091 

0.184±0.019 

3.12±0.16 

7.01±1.03 

-0.25 

12.39±0.60 

<0.9 

<-0.94 

J152116.59-I-001755 


0.039±0.0046 

0.274±0.014 

5.41±0.15 

9.51±0.70 

-0.60 

37.89±1.20 

1.19±0.28 

-0.90 

J154141.64-114409 


0.032±0.0077 

0.155±0.017 

2.91±0.16 

10.74±1.14 

-0.51 

34.52±1.13 

1.2±0.3 

-1.11 

J163426.87-172139 


0.039±0.0094 

O.lOlibO.018 

1.70±0.17 

3.57±1.15 

-0.42 

9.49±0.55 

<0.84 

<-0.63 

J164107.22-054827 


0.086±0.0083 

0.423±0.020 

3.14±0.15 

6.26±0.89 

-0.02 

6.62±0.48 

2.3±0.29 

-0.43 

J165305.40-010230 


0.083±0.0074 

0.191±0.015 

2.56±0.14 

5.31±0.93 

-0.36 

12.21±0.56 

<0.78 

<-0.83 

J165742.88-174049 


0.073±0.0102 

0.186±0.026 

2.82±0.24 

8.60±1.01 

-0.31 

17.48±0.71 

<0.78 

<-1.04 

J170204.65-081108 


0.021±0.0690 

0.074±0.053 

3.05±0.26 

12.32±1.40 

-0.74 

67.59±2.07 

<1.02 

<-1.08 

J170325.05-051742 


0.021±0.0082 

0.199±0.018 

2.35±0.24 

11.66±1.42 

-0.39 

28.77±0.96 

1.02±0.27 

-1.05 

J170746.08-093916 


0.119±0.0073 

0.342±0.020 

3.46±0.28 

3.27±1.26 

-0.52 

10.86±0.52 

<1.02 

<-0.51 

J193622.58-335420 


0.031±0.0069 

0.127±0.016 

2.34±0.14 

5.27±0.96 

0.00 

5.27±0.51 

1.86±0.36 

-0.45 

J195141.22-042024 


0.030±0.0178 

0.065±0.036 

2.55±0.15 

8.56±1.02 

-0.38 

20.52±1.09 

<1.05 

<-0.91 

J195801.72-074609 


0.056±0.0086 

0.203±0.018 

3.29±0.16 

7.44±1.06 

-0.64 

32.79±1.06 

<0.93 

<-0.90 

J200048.58-280251 


0.027±0.0169 

0.113±0.017 

3.21±0.17 

7.19±1.20 

-0.33 

15.33±0.66 

<0.96 

<-0.87 

J202148.06-261159 


<0.015 

<0.065 

1.03±0.15 

6.27±1.01 

-0.04 

6.82±0.55 

4.4±0.38 

-0.15 

J204049.51-390400 


0.070±0.0077 

0.254±0.017 

2.75±0.15 

4.02±0.91 

-0.44 

10.95±0.57 

5.1±0.43 

0.10 

J205946.93-354134 


0.052±0.0069 

0.182±0.015 

2.94±0.14 

4.75±0.99 

-0.28 

9.13±1.07 

<0.99 

<-0.40 


^ SDSS r-band AB magnitude 


redshifts for the sample to construct an “ensemble” rest- 
frame SED for our sample in Figure]^ We show ensemble 
SEDs of several comparison samples in Figure In all 
the plots we also show the tracks for several templates of 
nearby well-studied sources. The ensemble SEDs must be 
interpreted carefully because the choice of normalization 
wavelength affects the relative appearance of dispersion 
between the points at other wavelengths. We have cho¬ 
sen to normalize the templates and data at rest frame 4.6 
/rm. The selection of 4.6^m has the disadvantage that 
this spectral region may suffer significant extinction from 
a thick torus or other nuclear dusty structure, however 
it is the longest wavelength (and hence has the lowest 
optical depth) which both avoids the PAH and silicate 
features and which still lies within the WISE rest-frame 


wavelength range for all of the ALMA sources. 

The source t empla tes in all five figures are from 


Polletta et al. (2007), and include the cool starburst- 
dominated ULIRG Arp 220, the starburst M82, and the 
broad-line dusty QSO Mrk 231. The torus mod el in 
Figures [4||^ is based on the taper ed disk models of |Efs-| 
tathiou S Rowan- Robinsonl (1 19951) and has been htted to 


one o f our best sampled SEDs in this paper (see Section 


5.2.1). This torus model has an opening angle of 45° 


an inclination angle of 54° and a UV equatorial optical 
depth T,y=500. The intrinsic AGN SEDs, from |Shang| 


et al. (2011) and Mullaney et al. (2011), are empirical 
SEDs of nearby AgIN from which the host galaxy light 
has been subtracted, and they are quite similar to each 


other and to other published intrinsic AGN SEDs (Elvis 
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Table 2 

Herschel, JCMT and CSO Photometry 


WISE Name 

ALMA 

Redshift 

fro 

firo 

/250 

/sso 

fsoo 

/sso 


subsample 


PACS 

PACS 

SPIRE 

(mjy) 

SPIRE/CSO 

SPIRE 

JCMT/ALMA 

W0524+3005 



92.3±2.1 

119.8±7.3 

38.9±6.2 

<15.9 

<19.8 


W0526+1259 



31.3±2.0 

54.5±7.3 

34.3±6.3 

18.9±5.2 

<19.5 


W0537+3947 



33.7±1.7 

<26.1 





W0541+1130 



23.6±1.8 

<30.0 





W0844+7420 



106.1±2.0 

130.8±6.9 

62.0±6.5 

44.1±5.4 

24.6±6.2 


WlOOl-2141 






<45= 



W1025+6128 

W1331-3913 



<5.7 

<24.0 


30±13= 



W1332+7907 

W1343-1136 

ALMA 

2.49=" 

56.9±1.5 

59.6±7.8 


<45= 


2.34±0.31‘* 

W1400-2919 

ALMA 

1.67=" 




<45= 


<0.9=* 

W1500-0649 

ALMA 

1.500=" 

91.0 ±3.3 

171.6±7.4 




6.11±0.28‘* 

W1501+1324 


0.505*" 

199.9±1.7 

200.2±6.9 




< 6 . 6 = 

W1501+3341 



8.1±1.9 

49.4±7.0 





W1505+0219 



12 . 2 ± 1.6 

50.0 ±7.9 





W1517+3523 


1.515*" 

53.5±1.6 

69.8±8.5 




<5.7= 

W1921+7349 



19.4±2.0 

49.5±7.8 





W2005+0215 



14.2 ±1.9 

< 21.6 

<18.6 

<16.2 

<19.8 


W2059-3541 

ALMA 

2.380=" 

11.7±2.7 

<19.2 

<19.2 

<16.2 

<19.5 

<0.99=* 


BOO Table [3 
^ Palomar 200 inch 


CSO 

^ ALMA STOhm flux 
^ Jones et al. (2015) 


et al.||1994l 

Richards et al.| 

2006| INetzer et al. 

o 

o 

Aiq 

HD 

GG 

sef et al.|20^ 

ot a ~ 4 pc 

LO). In Figure [3 

radius dust toi 

we include the modeled fc 
us in the nearby AGN IN 


3081 that Ramos-Almeida et al. (2011) derived based 
on Herschel data and subarcsecond MIR imaging with 
Gemini T-ReCS. 

Most of our WISE-NVSS-ALMA sample is much more 
strongly MIR-dominated than the starburst templates 
(M82 and Arp 220) in Figure and also con ipared to 
the 1 < z < 3 Submillimeter Galaxies (SMGs) (Magnelli 


et al. 2012), which we illustrate to provide a compari¬ 
son to high redshift starburst-dominated systems. The 
WISE radio-blind Hot DOG samples (EWB12; Jo nes et 


al. 2014), a nd our northern sample of RP sources ( Jones, 
et al.|[Ml5l are similar to the ALMA sample in showing 


a low submm/MIR flux density ratio compared to the 
galaxy templates. Generally speaking, the WISE sources 
lie between the intrinsic SEDs and the Mkn 231 template, 
though some of the upper limits lie below even the intrin¬ 
sic SEDs shown. We also show some representative red- 


shift ^1-3 radio galaxies (Arch ibald et al. 2001 Grimes 
et al.|2005[|Seymo ur et al.|2007|) and broad-lin e (optically 
selected) quasars ( Priddey et al.||2003 2007) with avail¬ 
able 850/rm data ih the plot, however a detailed compar¬ 
ison to the FIR-submillimeter properties of well selected 
samples of AGN of various types is beyond the scope of 
this paper. The broad-line QSOs shown tend to have 
similar colors to our sample, while the RGs shown tend 
to be a bit more FIR-strong relative to the MIR, however 
these trends may be dominated by selection effects and 
should be used only as a very general comparison to our 
sample. 

The conclusions from Figure are emphasized in the 
ensemble SED (Figurel^top), where we see a strong simi¬ 
larity amongst the MIiGsubmm SEDs of the ALMA sam¬ 
ple. They resemble the high optical depth torus model 



Figure 3. Ratio of ALMA 870 //m to WISE 22 fim flux 
density vs. redshift for the ALMA sample (red points and 
limits). Histogram: sources without redshift (submillimeter- 
detected sources: filled bars; upper limits: empty bars). The 
comparison templates and source samples are described in the 
text. The three intrinsic AGN curves (purple dash-dot Ines) 
are the mean [log(L 2 —lOkev/erg s~^) < 42.9], meanfall], mean 
[log(L 2 -i 0 kev/ 6 rg s“^) > 42.9] SEDs of Mullaney et al'j (2011 \ 
in descending order which cover 6-1000 /im rest frame; the full 
dispersion of intrinsic AGN SEDs reported in that paper is signifi¬ 
cantly wider. The three ALMA sources lying on or above the Mkn 
231 template are W0417-2816, W0652-2006 and W0714-3635. *: 
The comparison samples include data at 850 fim and at 24 {im. 


in the rest frame NIR-MIR, as expected since this model 
was designed to fit one of our sources, however they are 
systematically steeper than any of the galaxy or intrinsic 
AGN templates in this wavelength range, when normal¬ 
ized at 4.6 fim. There is a lack of observational data 
available in the rest frame FIR, although the IRAS lim- 
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Table 3 

Redshifts, Luminosities and Radio Power 


WISE Name 

Redshift log Lagn 

(L©) 

Min Max 

30K 

log Lbb 

(Lo) 

50K 90K 

120K 

Min 

Og Lxotal 

(L©) 

Max 

Best 

Tduat 

(K)- 

Best model 

log PsGHz 
(WHz-i) 

Q22 

k-corr’d 

W0304-3108 

1.54 

13.01 

13.16 

11.49 

12.46 

12.92 


13.02 

13.36 

13.36 

90 

26.58 

-1.30 

W0306-3353 

0.78 

12.19 

12.23 

10.98 

12.07 

12.28 


12.25 

12.54 

12.53 

<50,90> 

24.92 

-0.55 

W0354-3308 

1.37 

12.70 

12.72 

<11.14 

<12.14 

<12.59 


<12.71 

<12.96 

<12.96 

90 

25.61 

-1.15 

W0404-2436 

1.26 

12.58 

12.63 

11.46 

12.48 

12.93 


12.61 

13.11 

13.11 

90 

25.67 

-1.69 

W0409-1837 

0.67 

12.73 

12.74 


<11.99 

<12.4 


<12.80 

<12.90 

<12.90 

90 

25.65 

-0.28 

W0417-2816 

0.94 

11.95 

12.01 

11.48 

12.55 

12.98 


12.12 

13.02 

12.47 

<30,50> 

25.54 

-1.38 

W0439-3159 

2.82 

13.44 

13.45 

12.04 

12.84 

13.33 

14.17 

13.47 

14.24 

13.70 

90 

26.76 

-1.12 

W0519-0813 

2.05'^ 

13.01 

13.02 

<11.30 

<12.2 

<12.68 

<13.58 

<13.03 

<13.68 

<13.18 

90 

26.53 

-1.09 

W0525-3614 

1.69 

12.51 

12.60 

<11.29 

<12.25 

<12.71 

<13.75 

<12.62 

<13.77 

<12.94 

90 

25.64 

-0.55 

W0526-3225 

1.98 

13.44 

13.54 

12.40 

13.31 

13.79 


13.57 

13.95 

13.95 

90 

27.33 

-1.18 

W0536-2703 

1.79 

13.05 

13.05 



12.94 


13.00 

13.30 

13.30 

90 

25.90 

-0.54 

W0549-3739 

1.71 

12.56 

12.56 

11.39 

12.34 

12.81 


12.59 

13.00 

13.00 

90 

25.97 

-0.92 

W0612-0622 

0.47 

12.03 

12.04 

10.89 

12.05 

12.20 


12.07 

12.43 

12.43 

90 

25.14 

-0.32 

W0613-3407 

2.18 

13.20 

13.20 

<11.43 

<12.31 

<12.79 

<13.68 

<13.21 

<13.80 

<13.34 

90 

26.57 

-0.90 

W0614-0936 

b 

12.94 

12.95 

<11.40 

<12.31 

<12.78 

<13.69 

<12.96 

<13.76 

<13.17 

90 

25.87 

-0.38 

W0630-2120 

1.44 

12.52 

12.57 

11.72 

12.71 

13.17 


12.63 

13.26 

12.82 

<30,50> 

25.86 

-1.33 

W0642-2728 

1.34 

12.36 

12.36 

11.34 

12.35 



12.40 

12.66 

12.66 

50 

25.50 

-1.42 

W0652-2006 

0.60 

11.58 

11.59 

11.11 

12.24 



11.71 

12.33 

12.12 

<30,50> 

24.74 

-0.39 

W0702-2808 

0.94'" 

12.04 

12.04 

<11.10 

<12.17 



<12.09 

<12.41 

<12.41 

50 

25.23 

-1.00 

W0714-3635 

0.88 

11.88 

11.89 

11.18 

12.26 



11.97 

12.41 

12.41 

50 

25.34 

-1.09 

W0719-3349 

1.63 

12.57 

12.57 


12.75 

13.22 


12.97 

13.31 

13.01 

<50,90> 

26.71 

-1.33 

W0811-2225 

1.11 

12.63 

12.63 



<12.65 


<12.94 

<12.94 

<12.94 

90 

25.75 

-1.53 

W0823-0624 

1.75 

13.11 

13.22 

<11.36 

<12.3 

<12.76 


<13.12 

<13.35 

<13.35 

90 

26.70 

-1.10 

W1308-3447 

1.65 

12.99 

12.99 

11.22 

12.18 

12.65 


13.00 

13.15 

13.15 

90 

26.74 

-1.40 

W1343-1136 

2.49 

13.00 

13.02 


12.43 

12.92 

13.78 

13.12 

13.85 

13.27 

90 

26.23 

-0.75 

W1400-2919 

1.67 

13.10 

13.11 

<10.78 

<11.73 


<13.16 

<13.11 

<13.43 

<13.43 

120 

26.63 

-1.13 

W1412-2020 

1.82 

13.02 

13.17 

11.52 

12.45 

12.92 


13.03 

13.36 

13.36 

90 

25.95 

-0.43 

W1434-0235 

1.92 

12.89 

12.89 

<11.08 

<12.00 

<12.47 

<13.39 

<12.90 

<13.51 

<13.03 

90 

26.62 

-1.22 

W1439-3725 

1.19 

12.29 

12.29 

<10.79 

<11.81 

<12.26 


<12.30 

<12.58 

<12.58 

90 

25.58 

-1.70 

W1500-0649 

1.50 

13.07 

13.50 


12.84 

13.33 


13.52 

13.59 

13.52 

90 

26.10 

-0.89 

W1510-2203 

0.95 

12.60 

12.61 

<10.74 

<11.81 

<12.23 


<12.62 

<12.76 

<12.76 

90 

25.57 

-0.84 

W1513-2210 

2.20 

13.26 

13.26 

11.86 

12.74 

13.22 


13.28 

13.54 

13.54 

90 

26.68 

-0.89 

W1514-3411 

1.09 

12.44 

12.48 

<10.81 

<11.86 

<12.29 


<12.45 

<12.70 

<12.70 

90 

25.57 

-1.36 

W1521+0017 

2.63'^ 

13.61 

13.61 



12.63 


13.65 

13.65 

13.65 

90 

25.60 

-1.04 

W1541-1144 

1.58 

12.81 

12.94 

11.13 

12.10 

12.57 


12.95 

13.01 

13.01 

90 

26.40 

-1.15 

W1634-1721 

2.08 

12.83 

12.83 

<11.08 

<11.98 

<12.46 

<13.36 

<12.84 

<13.47 

<12.98 

90 

26.11 

-0.80 

W1641-0548 

1.84 

12.94 

13.09 


12.41 

12.88 


13.05 

13.30 

13.30 

90 

25.83 

-0.37 

W1653-0102 

2.02 

13.00 

13.00 

<11.04 

<11.94 

<12.41 

<13.31 

<13.00 

<13.48 

<13.10 

90 

26.19 

-0.73 

W1657-1740 

b 

13.16 

13.17 

<11.03 

<11.94 

<12.42 

<13.32 

<13.17 

<13.55 

<13.60 

<90,120> 

26.34 

-0.68 

W1702-0811 

2.85 

13.60 

13.60 

<11.28 

<12.07 

_ ... 

<13.4 

<13.60 

<13.81 

<13.81 

120 

27.26 

-1.20 

W1703-0517 

1.80'^ 

13.11 

13.51 

11.12 

12.05 

12.48 

13.45 

13.12 

13.61 

3.60 

<90,120> 

26.91 

-0.74 

W1707-0939 

b 

12.98 

13.04 

<11.15 

<12.06 

<12.54 

<13.43 

<12.99 

<13.56 

<13.56 

120 

26.14 

-0.89 

W1936-3354 

2.24'" 

13.05 

13.18 

11.45 

12.32 


13.69 

13.06 

13.78 

13.78 

120 

25.64 

-0.40 

W1951-0420 

1.58 

12.80 

12.83 

<11.08 

<12.05 


<13.45 

<12.84 

<13.54 

<13.54 

120 

26.18 

-1.03 

W1958-0746 

1.80 

12.98 

12.98 

<11.08 

<12.01 

<12.48 

<13.41 

<12.99 

<13.54 

<13.10 

90 

26.51 

-0.99 

W2000-2802 

2.28 

12.75 

12.90 

<10.72 

<11.59 

<12.07 

<12.96 

<12.75 

<13.17 

<12.96 

90 

25.99 

-0.73 

W2021-2611 

2.44 

12.97 

13.51 

11.85 

12.68 


14.04 

13.00 

14.15 

14.15 

120 

25.93 

-0.45 

W2040-3904 

b 

12.92 

13.15 

11.85 

12.76 

13.23 


12.96 

13.49 

13.49 

90 

26.14 

-0.81 

W2059-3541 

2.38 

13.15 

13.37 



<12.46 

<12.87 

<13.15 

<13.42 

<13.42 

90 

26.26 

-0.69 


<T1,T2> denotes an average between these two model fits. 
z — l assumed if no spectroscopic redshift exists. 

Uncertain redshift: single line, blended lines or weak lines. 


its (red limits at 100 < Arest < 10 ^m) help to constrain 
the flux of many of our WISE-NVSS sample, ruling out 
SEDs like Arp 220 and M82. 

The three sources with highest submillimeteUMlR flux 
ratios (relative to the templates) from Figure!^ W0417- 
2816, W0652-2006 and W0714-3635 (yellow points in 
Figure]^ are the strongest candidates for possessing sig¬ 
nificant ongoing star formation. Their presence in our 
sample may be the result of a selection effect due to 
strong 6.7-7.7 /rm PAH features falling in the W3 fil¬ 
ter, or an 11.7 /im PAH feature falling in the W4 hlter, 
as can be seen in Figure]^ 

In Figure [^bottom we add the 10 sources with known 


redshift from Jones et al.|(2015|) (orange triangles), which 
are drawn from our northern WISE-NVSS sample. They 
all have upper limits at 850 fim from JCMT and they 
show very similar SEDs to the ALMA subsample. We 
also compare our radio-selected samples to the radio¬ 


blind WISE Hot DOG samples of EWB12 and of Jones 


et al. (2014). We see a very close similarity between the 


radio-selected (red and orange symbols) and radio-blind 
samples (green, blue, cyan and yellow filled symbols). 
The radio-blind Hot DOG samples have systematically 
larger redshifts than our sample, displacing the two sets 
of SEDs from each other in rest wavelength somewhat, 
however the radio-selected and radio-blind samples fall 
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Figure 4. Rest-frame ensemble SEDs for WISE-selected ultra red syste ms. Top; the 45 sour ces in our ALMA sample with known redshift, 
normalize d at 4 .6 //m, compared with the temp lates from the libra ry of|Polletta et aL|(|2007|l, the torus model fitted to one of our sources 
in Section|5.2.1| and the intrinsic AGN SEDs ofjShang et aL|(|2011|| (labeTTedTTQR^SOl QSOl). The colored triangles in the radio 

correspond to the matched color template. The three reddest objects, relatiye to Mkn 231, from Figure]^ are highlighted in yellow. The 
upper limits in the rest ~ 12-100 fim range are largely from IRAS, with a few from Herschel (see text for details). Seyeral sources lie at 
redshifts where the 22 fim band coincides with the 9.7 fim silicate absorption feature, which may partially account for the turnoyer in their 
spectral shapes. Note the intermediate radio power compared to classical (eyoly ed) RQ and RL Q SOs. when nor malized to 4.6 nm power . 
Bottom: Rest-fra me SEDs for the W ISE-selected red Hot DOG sources fro m IWu et al.| ||2012|| (green circles), |Eisenhardt et al.| (|2012|| 
(yellow diamonds), [Jones et al.| ( |2014[ | (cyan diamonds) and the WLABs from [Bridge et S! h pDTBl l (blue circles) com pared to ouF^ampl^ 
(red circles; orange triangles (Jones et al. 2015). Also shown is the reddest of rUe dust-reddened type 1 quasars from [Banerji etltT j polil l 
(cyan asterisks). 


within a continuous band. The different redshift selec¬ 
tion function for the two samples may be unrelated to 
radio power, instead being due at least in part to the 
redder {W2 — Wi) threshold of EWB12. This proba¬ 
bly eliminates sources with silicate absorption that falls 
into the 22/rm filter, and therefore favors sources with 


z >1.5-2. Jones et al. (2015) have suggested that the 


northern sources from our IN VSS-WISE sample that they 
observed at 850 /xm at the JCMT may show slightly less 
steep rest-frame SEDs than the EWB12 samples. This 
is an interesting possibility that requires further study. 

The WISE samples of EWB12 have more data from 
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Lonsdale et al. 



Figure 5. Rest-frame SEDs for several compa rison Svitzer and Herschel samples. Top: Rest-frame SEDs f or 152 Spitzer DOGs which 
were selected in a similar manner to our sample l|Polletta et al.|2008 |Melbourne et al.|2012||Sajina et al.|2012||. The extreme redness in the 
UV-optical for many of the Sajina et al. (20121 sources is due to~meir sele ction criteria, jjott om: Rest-trame SEDs for 61 high-redshift 
MIR-selected starburst-dominate d ULIRGs from |Lonsdale et al.|(|2009l and|Fiolet et al. ||2009|| (red symbols) . 16 Spitzer 70 /rm-selected 
galaxies from|Farrah et ah]||2007ll (blue symbols) and bl blVltls witn nerscnet data Irom Magnelli et al.|||2012[| (green symbols). 


Herschel in the rest MIR-FIR than we do for our sample, 
and their ensemble FIR SEDs tend to fall in the region 
occupied by the intrinsic SED^the torus template and 
the Mkn 231 template (Figure]^. The two WISE-NVSS 
sources for which we have Herschel data fall in the same 
region. All of these sources lie well below both M82 and 
Arp 220 in this wavelength region, emphasizing the result 
from Figure The limited Herschel data which we 
do have for our sample, and the IRAS limits, support 
a picture in which our WISE-NVSS sample have similar 


FIR-submm SEDs to the radio-blind WISE samples. The 


one Hot DOG with published radio data (Eisenhardt et 


al. 20121 has lower radio power than our sample, when 
normalized to 4.6 /rm power. 

It is interesting to compare the WISE-NVSS o bjects to 
the reddest known broad-line type 1 quasar from|Banerji| 


et al. (2014), ULAS J1234-I-0907 (z = 2.50), shown as the 
cyan asterisks in Figurefflbottom. The SED shape is very 
similar to the Hot DOGs through the MIR-submm but 
is much less red than any of the WISE sample in the rest 





























Jet Feedback in Obscured Radio-Loud Quasars 


11 


NIR 1-5 /im region. Unlike the WISE red sample, ULAS 
J1234-I-0907 follows the Mkn 231 template through ^ 1 
/rm, only droppi ng steeply at shorter wavelengths. The 


best fit model by Banerji et al.|(2014| has Ay = 6. This 
result emphasizes the litelihood that the steep red WISE- 
optical SEDs for the WISE-selected samples are caused 
by heavy obscuration. 


selected “power-law” DOGs dPolletta et al. 

2008 Mel- 

bourne et al.||2012| |Sajina et al.|2012), incfud 

mg the two 


(2008) 


samples 


ignif- 


icantly redder ensemble rest NIR slopes than most the 
Spitzer DOG samples. The WISE sources also appear to 
turn over into the EIR/submm at systematically shorter 
wavelength than the DOGs, as can be seen by comparing 
Figuresjfflbotto m and[5|-top, and as prev iously discussed 


by EWBl2 and Jones et al. (2014 2015). To first order, 
this is likely to indicate a larger ratio of star formation 
to AGN accretion in the Spitzer DOGs than in the Hot 
DOGs. 

Lastly we illustrate in Figure [^bottom the ensembles 
SEDs of several samples of Spitzer-selected starburst- 
dominated Ultra Luminous InfraRed Galax i es (ULIRGs 


1 Farrah et al. 

2007 

Lonsdale et al. 

2009 

Fiolet et al. 

2009), and SN 

IGs observed with Herscfiel (|Magnelli et 


Sion than either the WISE sources or the Spitzer-selected 
DOGs, relative to MIR emission, and bluer optical-MIR 
SEDs. 

Turning to the cm-radio emission there is a steady de¬ 
crease in the average radio/4.6 pm flux density ratio from 
our sample through the Spitzer DOGs to the SMGs / 
starbursts, although there is a lot of overlap between the 
latter two samples. This in entirely consistent with the 
selection in favor of bright radio sources in our sample, 
and the dominance of star formation in the SMGs and 
starbursts. 

In summary, the WISE-selected ultra red samples 
have very similar SED shapes from the rest-frame NIR 
through the FIR, with no obvious difference between the 
radio-selected samples (this paper and Jones et al. 2015) 
and the radio-blind Hot DOGs (EWB12; Jones 14) ex¬ 
cept that the radio-blind samples have a larger mean 
redshift, which may be caused by the different selection 
functions. Together, these sources are redder than any 
other known source type in the NIR-MIR, and most of 
them turnover into the FIR at higher frequencies than 
the Spitzer DOGs, starbursts and SMGs. 

4.3. Synchrotron Contribution to the ALMA fluxes 

Before addressing the possible range of star forma¬ 
tion rates in these sources, we first consider the possi¬ 
bility that a fraction of the 345 GHz flux is due to syn¬ 
chrotron emission. Since we have selected compact radio 
sources in radiatively efficient AGN there is the possibil¬ 
ity that the radio emission is beamed and that some of 
our sources are blazers. We briefly address the possibil¬ 
ity that some of the 345 GHz emission arises from non- 
thermal synchrotron emission associated with the radio 
sources here, however we defer detailed discussion of this 
topic to the next paper in our series, in which we present 
the high resolution X-band (8-12 GHz) imaging from the 
VLA. 


We can make some preliminary conclusions from the 
measured spectral indices across the 8-12 GHz VLA X- 
band, which are more reliable than indices derived from 
non-contemporaneous and non-beam matched 1.4 GHz 
NVSS data and the much later X-band imaging. We find 
that the majority of the sample has steep spectral indices 
between 8-12 GHz; 42 of them have 0*2 < —0.8 (27 have 
indices steeper than -1.0), characteristic of optically thin 
synchrotron emission and potentially consistent with be¬ 
ing Gigahertz Peaked Sources (GPS; < Ikpc in size with 
a synchrotron peak ~ 1 GHz) or Compact Steep Spec¬ 
trum sources (CSS; < 20 kpc in siz e with a synchrotron 
peak below 1 GHz) (O’Dea 1998). Several sources are 
also resolved or multiple on scales of 1-10 kpc. For most 
of these steep spectrum sources the synchrotron contri¬ 
bution to the 345 GHz flux is likely to be <10%. 

The remaining seven sources have flat or inverted 8-12 
GHz spectral indices, however two of these show a steep¬ 
ening spectral index between 12 and 20 GHz (from our 
limited VLA K-band imaging) thus they are also likely 
to be dominated by optically thin synchrotron emission 
(W0526-3225 and W0823-0642). Two of the remain¬ 
ing 5 sources (W0642-2728 and W1434-0235) have an 
ALMA measurement (a detection and a limit, respec¬ 
tively) which is well in excess (by a factor of 5 or more) of 
the extrapolated flat radio SED. The other three sources 
(W0536-2703, W1412-2020 and W1634-1721) require a 
synchrotron peak beyond 12 GHz to avoid exceeding the 
ALMA flux density detection or limit, and these are also 
probably not blazers unless they are exceptionally vari¬ 
able. They are more likely to be High Erequency Peakers 
(small < 100 pc sources with synchrotron peaks above 4 
GHz; Dallacasa et al. 2000) and it is possible that their 
345 GHz flux has a significant contribution from optically 
thin synchrotron emission. None of the flat or inverted 
spectrum sources has a plausible SED that can explain 
the WISE data as synchrotron emission from a blazar. 

5. ISM MASS, SED FITS AND DERIVED PARAMETERS 

As noted in the previous section, the simplest descrip¬ 
tion for the NIR-EIR SED shapes of the majority of our 
sample is that they resemble the intrinsic shapes of lo- 


emission (Elvis et al. 

19941 

Richards et al.|2006[ Netzer 

et al. 

2007 

lAsset et a 

fmro 

IMulfaney et af.|2011[|Shang 

et al. 

2011 

). It we interpret the SEDs in this fashion then 


emission stems from re-radiation of accretion disk energy 
by a nuclear torus or other dusty structure. In this sec¬ 
tion we explore, as an alternative, the range of plausible 
contributions from star formation that may be permitted 
by reasonable models of the SEDs. 

5.1. Interstellar Medium Mass 

We derive the ISM masses assuming that 100% of the 
345 GHz flux is thermal dust emission. As noted in 
the previous section the thermal fluxes could be overes¬ 
timated by <10% for most of the sample due to possible 
contributions from synchrotron emission, and by up to 
100% for 3 sources. We do not correct the data for this 
contamination because it is highly uncertain at this point 
in time. 


[Scoville et al. (|2014|) derive the gaseous ISM mass 
(Mhi + tWHa) of populations of distant galaxies in the 
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Lonsdale et al. 


COSMOS field using ALMA 870/xm data as a measure 
of cool dust mass. They use a local sample of well 
studied galaxies to show that the observed ratio of 850 
/rm specific luminosity to ISM mass Lysso/im/AfisM = 
1±0.23x10^° ergs s“^ Hz“^ ATq^ for low redshift spi¬ 
rals, with a dispersio n of a factor o f abou t 5. Using this 
empirical calibration, Scoville et al. (2014) derive the fol¬ 
lowing relation for a hux density measured at observed 
frequency Vohs (valid for Arest > 250/rm on the Rayleigh- 
Jeans tail where the emission will be optically thin) to 
derive ISM masses from ALMA data for SMGs (their 
equation 12, valid for a dust temperature of 25K): 


f^ab, ^ n oo AfiSM 

mJy 1 O^°M 0 


( ^obs 


( GpcA 

\ i^850/im j 

/ Fo 

1 di ) 


Here Fuj/Lo corrects for the departure of the dust 
emission spectrum from the Rayleigh-Jeans tail as the 
redshift increases from 0 (and the rest frequency ap¬ 
proaches the peak of the spectrum) an d di is the lu¬ 


minosity distance. (Scoville et al. 2014) adopt 250 
as the minimum acceptable rest-frame wavelength for 
derivations of ISM mass, thus we expect our results 
to provide reasonable estimates for sources with z<2.5, 
which represents 94% of our sample, with only small in¬ 
accuracies for the 3 sources with 2.5<z< 2.85. Frj is a 
function of a ssumed dust tempera ture, which was taken 
to be 25K by|Scoville et al. (2014). The constant 0.83 in 
the equation is also proportional to T^. We derive ISM 
masses (Table columns 13 & 14), for dust tempera¬ 
tures of 30 K and 90 K, using the appropriate values for 
the constant, Fq, and Frj. The majority of the dust in 
these systems is unlikely to be as warm as 90 K, there¬ 
fore these values provide fairly strong lower limits to the 
ISM masses when we have an 870 ^m detection. We es¬ 
timate the overall uncertainty to be a factor of 5, based 
on the unknown dust temperature and the overall dis¬ 
persion for local starburst systems in Figure 1 of Scoville 
et al. (2014). 

The ISM masses for = 30 K range from 0.8 to 
56 xlO^°Af 0 , with a median for the 26 ALMA-detected 
sources of 5.9 xlO ^^M^. This may be c ompared to the 
values derived by Scoville et al. (|2014 ) for their mid¬ 
plus high-redshft sample galaxy stacks from the COS¬ 
MOS field, which cover a similar redshift range as our 
sources. Their “IR Bright” sample has a median ISM 
mass of 11.91±0.77 xIO^^AJq, about twice that of our 
sample. 


5.2. Spectral Energy Distribution Fits 

The SEDs of most of the red WISE-NVSS sources are 
dominated by warm dust in the MIR with a strong de¬ 
cline into the submm, and many of them display spec¬ 
troscopic evidence of an obscured, high excitation, AGN. 
The minimum AGN luminosity can be reasonably well 
determined from the warm dust emission which domi¬ 
nates the WISE data. The total 1-000/rm luminosity is 
less well constrained, however, for most sources, because 
we lack data in the rest-frame FIR. We can model the 
emission by assuming a dominant dust temperature, T^, 
for the cooler dust that peaks in the FIR wavelength 
range, however the luminosity of such a component de¬ 
pends on the 4‘^ power of T^, and is uncertain by 2-3 


dex without measurements at the FIR peak of the dust 
SED. 

Two sources have available Herschel data, W1500- 
0649 and W2059-3541, and for these we find that the 
SED shapes fall well below the SEDs of M82 and Arp 220 
longward of the MI^ resembling the radio blind-selected 
HotDOGs (Figure ffl. We have constructed radiative 
transfer (RT) models for W1500-0649 (z=1.50), which 
has the most detections (7 bands in total) that include 
good constraints on the peak of the SED, in order to pro¬ 
vide some insights into the nature of these sources. For 
the remaining sample there are insufficient data points to 
justify RT models, therefore we use a parametric “torus” 
model plus modified blackbody (BB) fits for the cooler 
dust. 

5.2.1. Radiative Transfer Model for WISE 1500-0649 

We have constructed example models for W1500-0649 
using both a tapered disk model (Efstathiou & Rowan- 
Robinson 1995) for the torus (see also Efstathiou et 
al. 2013) and the two-phase clumpy torus models of 
Stalevski et al (2012). For the cooler dust component 
radiating in the FIR we use the starburst models of Efs¬ 
tathiou, Rowan-Robinson & Siebenmorgen (2000) which 
were revised by Efstathiou & Siebenmorgen (2009). 

The tapered disk models have 4 parameters plus a nor¬ 
malization factor and the starburst models have 3 param¬ 
eters plus a different normalization factor. In the tapered 
disc AGN models we fix the opening angle of the torus 
©0 at 45° and vary the equatorial optical depth of the 
torus in the range 250-1250 {Ay « 50 — 250), the 
inclination i in the range 45 to 90° and the ratio of outer 
torus radius to inner radius rijrx from 20 to 160. In the 
starburst models we fix the initial optical depth of the 
molecular clouds Ty at 75 (which is the average value in 
the model grid) and the e-folding time of the starburst at 
20 Myr which from previous work appears to be a rea¬ 
sonable timescale for starbursts (e.g. Efstathiou et al. 
2000). In the starburst models we only vary the age of 
the starburst t*. 

We use a standard minimization technique to find 
the model parameters that best fit the data. =The best 
tapered disk fit is shown in Figure [^top (the AGN torus 
model is plotted with a blue dotted line and the starburst 
model with a red dashed line) assumes the following pa¬ 
rameters: = 500, rijrx = 160, i = 54° and U = 15 

Myr. 

The AGN torus has a derived luminosity of 1.9 x 
lOl^Lg. This needs to be multiplied by the anisotropy 
correction factor A (Efstathiou 2006) which for this par¬ 
ticular combination of parameters is 0.84 to give an AGN 
luminosity of 1.6 x IO^^Lq. The starburst luminosity is 
1.8 X IO^^Lq so the total luminosity of the system is 
predicted to be 3.4 x IO^^Lq. 

In Figure [^bottom we show the SED of W1500- 
0649 fitted with the two-phase clumpy torus models of 
Stalevski et al (2012) in combination with the starburst 
models of Efstathiou et al. (2009). We find that a good 
fit can be obtained with a torus that assumes a half¬ 
opening angle of 50 and an inclination of 90 degrees. The 
AGN luminosity is predicted to be 1.2 x IO^^Lq and the 
starburst luminosity 2.4 x lO^^L©. The starburst is still 
predicted to be a young system with an age of 10 Myr. 
We conclude that irrespective of the uncertain dust ge- 
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ometry, in W1500-0649 the AGN and the starburst emit 
comparable luminosity. 


5.2.2. Three-Component SED Fits 

We model the full sample by fitting the minimum num¬ 
ber of simple spectral components to the SEDs that will 
define a reasonable maximum star formation rate that is 
consistent with the observed data, and to provide min¬ 
imum and maximum estimates of the luminosity of the 
AGN-heated, dust. A third SED component in the opti¬ 
cal is used to estimate stellar mass, and we assume the 
AGN is completely obscured in the optical-UV. The re¬ 
sults are presented in Tables|3]^and illustrated in Figure 
(which is continued in the Appendix). A more thor¬ 
ough analysis of the SEDs using radiative transfer models 
for the full sample is in preparation by A. Efstathiou et 
al. 

The method fits a parametrized “torus” model to the 
MIR data and a modified blackbody (BB) to the longer 
wavelength SED. The “torus” model could describe a 
classical torus, or some other dusty structure heated by 
the AGN, including a spherical “cocoon” with 100% cov- 


J1 50048.73-064939.8 



Log (X//um) 


Figure 6. Best fit radiative transfer models for WISEJ150048.73- 
064939.8. Top: AGN tapered disk model of|Efstathiou & Rowan-| 
|Robinson| |jl995 \ which has = 500, r 2 lr\ = 160, i = 54“ and 
t, = 15~M vr.: Bottom: AGN Clumpy torus model of |Stalevski et| 
[ah] ||2012|| (blu e dotted lines’). The sta rburst models (red dashed 
lines) are from [Efstathiou et al.||2009|l. 


ering factor of warm dust, flwD- For the BB component 
the characteristic dust temperature is undefined for most 
sources, therefore we construct four models with different 
fixed dust temperatures. Modeling a range of dust tem¬ 
peratures within an individual source is not justified by 
the available data points. We do not expect all sources 
to be fitted well for each of these temperature choices, 
particularly the higher values, and we carry forward into 
the analysis only those SED fits that are viable. 

The methodology follows Sajina et al. (2012| but is 
simplified from their four dust components to only two, 
due to the limited SED data available for our sample. 
The AGN MIR emission is modeled with a parameteri- 
za tion which is consiste nt with the clumpy torus models 
of Nenkova et al. (20081. It has the functional form: 


u = 


(^)ae0.5. 


(JL)-0.5 + (JL)-3 

V 1^0 ' V I/Q ^ 


17W3 (2) 


where a and vq are free parameters. The MIR fit is 
phenomenological and so it is not characterized by spe¬ 
cific values for the dust temperature range, orientation, 
torus size or optical depth. We refer to this component 
henceforce as the “AGN” component. We emphasize 
that it could represent a structure of different shape than 
a classical torus, such as a more spherical cocoon, or a 
dusty NLR or polar wind. We also do not rule out in our 
later discussion that some of this warm emission could 
be contributed by a young compact starburst. 

The modified blackbody component has a fixed dust 
temperature and emissivity, /?. With only the single 
ALMA data point long ward of 22 /im for most sources, 
it is impossible to fit the long wavelength portion of the 
SED using free parameters, therefore a single greybody 
is the best approximation. This component may repre¬ 
sent dust heated by star formation, or dust heated by the 
AGN that is cooler than the “AGN component”, and of 
course would in reality it have a wide range of dust tem¬ 
peratures compared to the single value used here. The ht 
to the MIR region of the spectrum with the AGN compo¬ 
nent varies very little as the BB component temperature 
is changed, because the AGN component it is very well 
defined by the WISE data on the short wavelength side. 
For the four sources without measured redshift we as¬ 
sume z = 2. 

For the stellar luminosity and mas s we use either a 
100 or 600 Myr stellar population from Maraston (2005). 
This is constrained only by the i?-band data point in 
most cases, except for the two SDSS-detected sources, 
and the shape of the MIR SED at the shortest wave¬ 
lengths. 

We have made four fits to each source, each with a 
different temperature for the BB component. We have 
allowed the temperatures to take on a wide range so that 
we can interpret the luminosity as arising from disk-like 
distributed star formation (the coolest dust temperature, 
30 K); a starburst sim ilar to those found in local LIRGs 
and ULIRGs (50 K) ([Melbourne et al. 2012 Bendo et 
and two additional warmer temperatures that 
might be appropriate for very young and compa ct star- 
bursts and/ or for additional AGN-heated dust (Wilson 
et al.[[20l4| . For j3 we have selected a value of 1.5, con¬ 
sistent witn the range of values found in the literature. 
For the 50K model only, we instead used a value of 13=2, 
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consistent with the largest values found in the literature, 
in order to illustrate that the uncertainties due to the un¬ 
known dust temperature distribution in each source far 
outweigh the effect of the choice of /?. This may be seen 
by comparing the 30 K and 50 K model fits in Figure [8 
We do not expect any of these fits to be unique for any o' 
our sources; instead we use them to constrain the range 
of plausible AGN and starburst luminosities. 

To select appropriate dust temperatures for the two 
warmer BB models, we ran a set of models where was 
allowed to be a free parameter for the two sources which 
have Herschel data, and the results are shown in Figure 
fright, compared to the 30 K and 50 K fixed-T fits in 
^e left panels. For W1500-0649, which was fitted in the 
previous section with radiative transfer models, there are 
two plausible fits. The 50 K BB model is reasonable, al¬ 
though the IRAS 60 ^m limit is slightly exceeded. This 
fit requires a large fraction of the FIR emission to be ex¬ 
plained by the warm AGN component (dashed cyan line). 
The best fit free-T^ model has = 89 K and provides a 
total MIR-submm luminosity of log(LTotai/L0) = 13.52, 
and a division between the AGN and BB components of 
1.17 and 2.14 xIO^^Lq. This is very similar result to the 
radiative transfer model. For W2059-3541 a dust tem¬ 
perature greater than 50 K is required to match the 70 
fim PAGS data. The best-fit T^-free model has Td = 120 
K, with a total MIR-submm luminosity of log Lboi/L© < 
13.33 (the 870/im measurement is a non-detection for this 
source), assuming /3=1.5. 

For the remainder of the sample, for which free-T^ fits 
are insufficiently constrained, we have adopted the two 
dust temperature values which have fitted W1500-0649 
and W2059-3541 successfully: 90 K (rounding up from 89 
K) and 120 K, for the two warmest of the 4 BB models. 
The two fits with the cooler dust temperatures, 30 K and 
50 K, are shown in the left hand panels of Figurej^ while 
the fits with the two warmer temperatures, 90 K and 120 
K, are shown in the right hand panels. The results are 
tabulated in Table [Sl Since the AGN luminosities do 
not vary much between the four models, we list only the 
minimum and maximum AGN component luminosities in 
columns 3 and 4 of Table s -Z^AGN-Min and LAGN-Max- 

Many of the BB fits are not successful in fitting all the 
data points; in particular the 120 K fits exceed the IRAS 
limits for many sources. Also the stellar component is not 
always fitted well by the code. We have not attempted 
to refine the fitting procedures to improve this situation 
because the current optical photometry is inadequate for 
this purpose. Only fits judged to be acceptable are tab¬ 
ulated and used in the subsequent analysis. 

For the two warmer BB models, the IRAS 60 /rm 
limit (and occasionally the 100 /rm limit) is helpful in 
constraining the dust temperature of the fits, for many 
sources. Four sources fail to find a reasonable fit for the 
IRAS limits with Td = 90 K, while for 29 sources, the 
IRAS 60 fjim limit rules out the 120 K model. In one case, 
W1500-0649, the 120 K model is ruled out not by IRAS 
but by ALMA and Herschel. In some cases the low tem¬ 
perature BB models fail to allow the overall model to fit 
the WISE data; these sources are better fit with one or 
both of the two warmer models. 

We derive star formation rates from the BB compo¬ 
nents of the models using the Kennicutt (1998) conver¬ 


sion from far-infrared luminosity to star formation rate: 

SFR = 4.4 10-3^Lbb(W) M© yr"^ (3) 

The SFR results are given in Table columns 3-6, for 
those models which achieved successful fits. 

We also list in Table [3] columns 9 and 10 the min¬ 
imum and maximum summed luminosity from the up 
to 4 viable LTotai (= Lagn + Lbb) models, LTotai-Min 
and Axotai-Max, for each source. In addition we list 
LTotai-Best in column 11, which represents the Axotai fit 
that best resembles the SED shapes of the two sources 
with well fitted SEDs: W1500-0649 and W2059-3541, 
and the Hot DOGs with well-sampled SEDs (Figure]^. 
The best fit model is indicated in column 12. In some 
cases the best match is derived from the average of two 
of the Axotai models. 

The luminosities derived for the AGN from the AGN 
model component range from log(LAGN-Min/.^0) = 
11.58 to log(LAGN-Max/L0) = 13.61. We also consider 
a maximum AGN luminosity derived from the Total- 
Best models: log(LTotai-Best/L©) = 12.01-14.15. The 
BB component luminosity depends strongly on assumed 
dust temperature. Gonsidering all acceptable fits to the 
SEDs for all of the 870 /rm-detected sources, the accept¬ 
able range in log^L^B/Le) is 10.89-14.17, and the corre¬ 
sponding SFR range is 13.5-25700 M© yr“^. The total 
summed AGN-I-BB luminosity range is log(LTotai/A0) = 
11.71-14.24. These results are summarized in Table SI 

5.3. Stellar Masses 

The stellar population fit can be used to constrain the 
host galaxy mass from the rest frame H-band absolute 
magnitude of the fitted stellar component. The host 
galaxy mass has only a small dependence on the stel¬ 
lar population selected by the fitting code. The masses 
do not depend on the dust temperature assumed for the 
BB dust component because the optical fit is dominated 
by the AGN component fitted to the WISE data. In 
a few cases, a self-consistent fit could not be obtained, 
indicating some possible confusion in the R-band data 
point. The H-band absolute magnitudes are upper lim¬ 
its in more than half the sample since there is no R- 
band measurement or detection, therefore most of the 
stellar mass estimates in Golumns 11-12 of Tableware 
also upper limits. The stellar masses derived from the 
models range from log(Mstars/Af0) = 9.83 — 11.09 for 
i?-band detected sources. The upper limits ranges from 
l0g(Mstars/A/©) < 9.82 to <11.34. 

The modeled intrinsic H-band absolute magnitudes 
and the derived stellar masses are highly uncertain and 
should be viewed as indicative only. We have no measure¬ 
ment of extinction for the optical or WISE data points, 
and no color information to help constrain the stellar 
populations. We also have little information on morphol¬ 
ogy, and therefore on the stellar components included in 
the flux density. It is also possible that some of the emis¬ 
sion is scattered AGN light. The fitting of the stellar 
component will be greatly improved by the use of the J 
and Ks data therefore we limit discussion of these results 
in this paper. 

5.4. Radio Powers, Radio Loudness and the q 
Parameter 
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Figure 7. SED models for the two sources for which we have sufficient data near the peak of the SED, from Herschel, to obtain well 
constrained fits across the MIR and FIR: W1500-0649 and W2059-3541. “AGN” (dashed cyan line); modified blackbody “BB” dust 
component (green lines); stellar population (red lines); summed model (black line). For the AGN and stellar components, and for the total 
fits, only one fit is shown in each panel, for clarity. The ALMA data point is the right-most point, and the limits near the peak of the cool 
component are 60 /im and 100 fim IRAS data. Left: models with = 50 K, = 2 (upper) and = 30 K, /3 = 1.5 (lower). Right: 
models in which the temperature of the BB dust component is allowed to float, for these two sources only, with 0 = 1.5. 


Table 4 

Minimum and Maximum Model Range in Luminosity, Mass, SFR &: Accretion Rate (M) 


Parameter 

AGN 

BB 

Total 

-870/im Detected 

Total-Best 

BB 

Total 

0/im Not Detected- 

Total-Best 





Min Value 

W0652-2006 

W0612-0622 

W0652-2006 

W0612-0622 

W2000-2802 

W0702-2808 

W0702-2808 

Max Value 

W1521-I-0017 

W0439-2159 

W0439-2159 

W2021-2611 

W0525-3225 

W1702-3225 

W1702-3225 

Redshift 

0.60-2.63 

0.47-2.82 

0.60-2.82 

0.47-2.28 

2.28-1.69 

0.94-2.85 

0.94 -2.85 

log L (L©) 

11.58-13.61 

10.89-14.17 

11.71-14.24 

12.07-14.15 

<10.72-<13.75 

<12.09-<12.46 

<12.41-<13.81 

SFR (M©/yr) 


13.5-25700 

22-25700 

a 

<9-<9800 

<22-<4360 

a 

log Mbh (M©)'^ 

7.66-9.69 


b 

8.11-10.23 


b 

<8.49-<9.89 

M (M©/yr)<^ 

0.24-26 


b 

0.7-90 


b 

<1.6-<41 


No star formation is present in the Total-Best mode! by definition, since it is designed to resemble the intrinsic AGN SEDs. 

Some Total=AGN-|-BB fits are inconsistent with a torus-like or intrinsic AGN SED shape; the Total-Best model is the preferable maximal fit for an AGN. 
A^bh accretion rate are directly proportional to the luminosity since we have assumed a fixed Eddington ratio and accretion efficiency. 


The rest frame 3 GHz radio powers in column 13 of Ta¬ 
ble^ are found to lie in the range log {P 3 Giiz/WHz~^) 
= 24.74-27.33, adopting a power law spectral index 
a = —1.0 {fi, oc z/“) for the /c-correction to 3 GHz 
rest frequency. The median/mean value is 25.97/26.05. 
We also list the rest-frame (i.e. fc-corrected) 922 = 
log(/22pm//i.4GHz) values in column 14. The MIR k- 
correction for q obviously depends strongly on the as¬ 
sumed rest frame SED. In our situation this is partic¬ 
ularly tricky because the unknown depth of the silicate 
feature will affect the observed /22 flux density strongly 
at redshifts near 1.5 (e.g. see Figure 1^. We have u sed 
the QS02 (“Torus”) template oflPolleto et al.|(|2007|) be¬ 
cause it was successfully fitted by them (their Figure 9) 
to a very red, Compton thick, Spitzer-selected obscured 
QSO that is similar in NIR-MIR spectral shape to our 
sources. To this we have grafted on the silicate absorp¬ 


tion fe^ure from the Arp 220 template of |Polletta et| 
al. (2007). We also derived the /c-correction without the 
added silicate feature. All of the templates have been 
convolved with the WISE 22 n m filter as a function of 
redshift by Polletta et al. (2007). Our quasars have com¬ 
parable fc-corrected ^22 values to the most radio-powerful 
and most radio-loud sources in the la rge Spitzer SXDF 
(Subarau X-ray Deep Field) sample of Ibar et al. (2008), 
and all of our sources would be considered to be mod- 
erately to very radio-loud after the fc-correction, based 
on this criterion. This is the case even if the large sil¬ 
icate optical depth is omitted from the template. This 
conclusion is also apparent from Figure [5 in which our 
quasars lie between the radio-quiet and radro-loud quasar 
templates, when normalized at 4.6 /im rest. 

5.5. Black Hole Masses and Accretion Rates 
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Figure 8. Example of the SED models for each source: data and 
fits as in Figur e[7| The full set of SED fits is shown in the appendix. 
Many of the BBmodels do not provide good fits, and these are not 
considered further in our analysis. In some cases no consistent fit 
is found for the stellar component. Left: models with = 50 K, 
/3 = 2 (upper) and = 30 K, j3 = 1.5 (lower). Right: models 
with 90 K (lower) and 120 K (upper), with 0 = 1.5; temperature 
choices based on the results for W1500-0649 and W2059-3541 in 
Figure [ t] ( continued in Appendix). 

Lacking a high quality spectroscopic indicator of BH 
masses such as the Mg II emission line width, we de¬ 
rive BH masses and accretion rates from the AGN lumi¬ 
nosities. These three parameters are therefore directly 
proportional to each other in this work. We assume an 
Eddingt on ratio, Af.hhi of 0.25, which is typical of z~2 
quasars (Kormendy fc Ho||2013[): 


AEdd = {4:TrGmpc/aT)MBH = 3.3 10"^ M^h (4) 

We assume a covering factor of unity and spherical 
symmetry, consistent with late stage mergers with heav¬ 
ily obscured nuclei. The covering factor cannot actually 
be this high in most sources, since emission lines from 
the NLR are visible in many cases; therefore our AGN 
luminosities and BH masses may be under-estimates due 
to missing emission in the X-ray to optical range which 
is not absorbed by the dust. We present the results in 
Tablecolumns 3-5, using our Lagn-miii, LAGN-Max, 
and LTotai-Best estimators. The median/mean BH mass 
values are 1.0/1.15, 1.07/1.45, and 1.55/2.75 xlO^M© 
respectivel y. In columns 6-7 we also list the BH masses 
derived by Kim et al. (2013), based on the [O III]A5007A 
line luminosity. 

The accretion rates are derived from the AGN lumi¬ 
nosities assuming an efficiency fo r the conversion of mat - 
ter into radiant energy of e = 0.1 (Heckman & Best 2014): 


Lboi = eMc2 (5) 

The rates are listed in columns 8-10 of Table and 
range from 0.24 to 25.3 A/gyr”^, with median/mean val¬ 
ues of 6.1/5.6, 7.6/6.2 and 14.5/9.1 MQyr“^ respectively 


for the AGN-Min, AGN-Max and Total-Best values. The 
overall ranges in BH mass and accretion rate are listed 
in Table m 

6. DISCUSSION 

Our goals are to search for the most luminous obscured 
quasars at redshifts ~l-3, the peak epoch of massive BH 
building, which are in the process of quenching star for¬ 
mation, and to investigate the role of radio jets in that 
process. We are specifically interested in the kinetic role 
of moderate to high power jets on the ISM within the 
galaxy host, while the AGN is still accreting strongly in 
“quasar-mode”. This is in contrast to the role of jets 
in typical “radio-mode” AGN, which are thought to be 
accreting at low rates and to have a role in maintaining 
galaxies free of infalling gas. As outlined in Section 
models which take the porosity of the ISM into account 
show that high power jets could be quite effective in a 
dense dusty environme nt, such as found the central re- 


gions of major mergers (Wagner & Bicknell 2011 

-— 


Wagner 


We have shown that WISE has found sources that are 
steeper (redder) in the rest-frame 1-10/rm range than 
most known Spitzer-selected red sources, including the 
so-called DOGs. The two other samples of very red 
WISE sources, the Hot DOGs and the WISE Lyman Al¬ 
pha Blobs (WLABs) (EWB12) which were selected with¬ 
out regard for radio emission, have extremely red rest- 
frame 1-10 /im SEDs which are very similar to those our 
our radio-selected sample. The radio-blind samples have 
higher average redshifts than our radio powerful sam¬ 
ple, probably as a result of the different WISE-optical 
color selection criteria. We attribute our ability to iden¬ 
tify these extremely red sources to the much larger sur¬ 
vey volume of WISE compared to Spitzer. These sys¬ 
tems are candidate luminous quasar-mode AGN which 
are highly obscured, and the radio-bright ones selected 
here potentially have young jet activity. Recent X-ray 
observations of a few Hot DOGs confirm the likelihood 


of highly buried AGN ( 

Stern et al. 

2014 

Pinconcelli et al. 

2U151. 'i'he hydrogen column tor V 

^1835+4355 was tound 


We also find that the WlSE-iN VSS-ALJVIA sample 
sources are more strongly dominated by AGN than the 
Spitzer DOGs. Star formation rates of hundreds up to 
a few thousand Mq yr“^ could also be present in some 
of the systems, although the IR-submm SEDs of some of 
the ALMA sources, in particular those with only upper 
limits at 345 GHz, may be consistent with the AGN torus 
model in Eigure without any SF contribution. There¬ 
fore these sources are indeed objects in which it is likely 
that the accretion rates are still very high but the star 
formation rates are low relative to the accretion power, 
and could thus be ideal sources for investigating recent 
and ongoing quenching by jet-powered AGN feedback. 

Some of our systems could be High Excitation Radio 
Galaxies (HERGs) seen at high inclination through an 
optically thick torus, and it is possible that we have in¬ 
cluded some with lobes that are unresolved by the 45" 
NVSS beam. We will show in the next paper (Garol 
Lonsdale et al. 2015, in preparation) that our 8-12 GHz 
VLA data rule out this scenario for most sources, al¬ 
though a small subset of the VLA sample of 156 red ob¬ 
scured quasars does indeed turn out to have large (several 
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Table 5 

Masses and Accretion Rates 


WISE Name 

Redshift 

Min 

log Mbh 
(Mq) 

Max Total- 
Best 

log MBH:Kim,ctoL 

(Mq) 

Lower Upper 

Accretion Rates Absolute log Mgtars 

(Mo/j/r) H- (Mq) 

Min Max Total- Magnitude 

Best 

log Mism 
(Mq) 

30K 90K 

W0304-3108 

1.54 

9.09 

9.24 

9.44 



6.52 

9.20 

14.50 

-24.5 

10.6 

10.83 

10.20 

W0306-3353 

0.78 

8.27 

8.31 

8.61 



0.99 

1.08 

2.18 

-23.9 

10.4 

10.51 

9.96 

W0354-3308 

1.37 

8.78 

8.80 

<9.04 



3.19 

3.34 

<5.82 

-23.6 

10.2 

<10.52 

<9.91 

W0404-2436 

1.26 

8.66 

8.71 

9.19 



2.42 

2.72 

8.14 

-25.1 

10.8 

10.86 

10.26 

W0409-1837 

0.67 

8.81 

8.82 

<8.98 



3.42 

3.50 

<5.10 

-22.6 

9.8 

<10.44 

<9.90 

W0417-2816 

0.94 

8.03 

8.09 

8.55 

9.06 

9.72 

0.57 

0.65 

1.87 

-22.9 

9.9 

10.63 

10.06 

W0439-3159 

2.82 

9.52 

9.53 

9.78 



17.54 

17.95 

31.56 

-23.6 

10.3 

11.15 

10.38 

W0519-0813 

2.05 

9.09 

9.10 

<9.26 



6.52 

6.67 

<9.72 

>-23.1 

<10.0 

<10.53 

<9.85 

W0525-3614 

1.69 

8.59 

8.68 

<9.02 



2.06 

2.54 

<5.58 

-23.5 

10.2 

<10.59 

<9.94 

W0526-3225 

1.98 

9.52 

9.62 

10.03 



17.54 

22.08 

56.80 

-24.1 

10.4 

11.64 

10.96 

W0536-2703 

1.79 

9.13 

9.13 

9.38 



7.14 

7.14 

12.69 

-23.7 

10.3 

10.82 

10.16 

W0549-3739 

1.71 

8.64 

8.64 

9.08 



2.31 

2.31 

6.42 

>-22.7 

<9.9 

10.69 

10.04 

W0612-0622 

0.47 

8.11 

8.12 

8.51 



0.68 

0.70 

1.71 

-23.8 

10.2 

10.53 

10.01 

W0613-3407 

2.18 

9.28 

9.28 

<9.42 

9.07 

9.46 

10.09 

10.09 

<14.02 

>-22.6 

<9.8 

<10.64 

<9.94 

W0614-0936 

a 

9.02 

9.03 

<9.25 



5.55 

5.68 

<9.51 

-22.9 

<9.9 

<10.64 

<9.96 

W0630-2120 

1.44 

8.60 

8.65 

8.90 



2.11 

2.37 

4.17 

-24.0 

10.4 

11.08 

10.46 

W0642-2728 

1.34 

8.44 

8.44 

8.74 



1.46 

1.46 

2.88 

-24.1 

10.4 

10.72 

10.11 

W0652-2006 

0.60 

7.66 

7.67 

8.20 

8.85 

9.61 

0.24 

0.25 

0.84 

>-21.3 

<9.2 

10.70 

10.16 

W0702-2808 

0.94 

8.12 

8.12 

<8.49 

8.92 

9.73 

0.70 

0.70 

<1.64 

-21.9 

9.5 

<10.58 

<10.01 

W0714-3635 

0.88 

7.96 

7.97 

8.49 

9.26 

9.34 

0.48 

0.49 

1.64 

-21.7 

9.5 

10.69 

10.13 

W0719-3349 

1.63 

8.65 

8.65 

9.09 



2.37 

2.37 

6.47 

-22.5 

9.8 

10.51 

9.74 

W0811-2225 

1.11 

8.71 

8.71 

<9.02 



2.72 

2.72 

<5.56 

-24.0 

10.4 

<10.61 

<10.02 

W0823-0624 

1.75 

9.19 

9.30 

<9.43 



8.20 

10.57 

<14.23 

-24.2 

10.5 

<10.64 

<9.99 

W1308-3447 

1.65 

9.07 

9.07 

9.23 



6.22 

6.22 

9.07 

-24.4 

10.6 

10.53 

9.89 

W1343-1136 

2.49 

9.08 

9.10 

9.35 



6.37 

6.67 

11.81 

-25.7 

11.1 

10.74 

10.00 

W1400-2919 

1.67 

9.18 

9.19 

<9.51 



8.02 

8.20 

<17.22 

-24.8 

10.7 

<10.09 

<9.44 

W1412-2020 

1.82 

9.10 

9.25 

9.44 



6.67 

9.42 

14.72 

>-24.7 

<10.7 

10.80 

10.14 

W1434-0235 

1.92 

8.97 

8.97 

<9.11 



4.94 

4.94 

<6.82 

-24.5 

10.6 

<10.34 

<9.67 

W1439-3725 

1.19 

8.37 

8.37 

<8.66 



1.24 

1.24 

<2.40 

>-23.2 

<10.1 

<10.21 

<9.61 

W1500-0649 

1.50 

9.15 

9.58 

9.60 



7.48 

20.14 

21.10 

>-23.2 

<10.1 

11.18 

10.56 

W1510-2203 

0.95 

8.68 

8.69 

<8.84 



2.54 

2.59 

<3.68 

>-23.1 

<10.0 

<10.23 

<9.66 

W1513-2210 

2.20 

9.34 

9.34 

9.62 



11.59 

11.59 

22.16 

>-23.9 

<10.4 

11.07 

10.36 

W1514-3411 

1.09 

8.52 

8.56 

<8.78 



1.75 

1.92 

<3.16 

>-23.6 

<10.2 

<10.25 

<9.67 

W1521+0017 

2.63 

9.69 

9.69 

9.73 



25.94 

25.94 

28.66 

>-26.1 

<11.2 

11.75 

11.20 

W1541-1144 

1.58 

8.89 

9.02 

9.09 



4.11 

5.55 

6.48 

>-23.5 

<10.2 

10.46 

9.83 

W1634-1721 

2.08 

8.91 

8.91 

<9.06 



4.31 

4.31 

<6.14 

>-25.1 

<10.8 

<10.31 

<9.62 

W1641-0548 

1.84 

9.02 

9.17 

9.38 



5.55 

7.83 

12.66 

>-26.4 

<11.3 

10.75 

10.09 

W1653-0102 

2.02 

9.08 

9.08 

<9.18 



6.37 

6.37 

<8.00 

>-25.4 

<11.0 

<10.28 

<9.60 

W1657-1740 

a 

9.24 

9.25 

<9.68 



9.20 

9.42 

<25.42 

>-25.1 

<10.8 

<10.28 

<9.60 

W1702-0811 

2.85 

9.68 

9.68 

<9.89 



25.35 

25.35 

<41.35 

>-25.2 

<10.9 

<10.38 

<9.61 

W1703-0517 

1.80 

9.19 

9.59 

9.68 



8.20 

20.61 

25.62 

>-25.1 

<10.8 

9.92 

9.16 

W1707-0939 

a 

9.06 

9.12 

<9.64 



6.08 

6.98 

<23.22 

>-25.1 

<10.8 

<10.39 

<9.71 

W1936-3354 

2.24 

9.13 

9.26 

9.86 



7.14 

9.64 

38.33 

>-24.5 

<10.6 

10.97 

10.32 

W1951-0420 

1.58 

8.88 

8.91 

<9.62 



4.02 

4.31 

<21.97 

>-24.9 

<10.7 

<10.40 

<9.76 

W1958-0746 

1.80 

9.06 

9.06 

<9.18 

9.53 


6.08 

6.08 

<8.00 

>-24.7 

<10.7 

<10.36 

<9.70 

W2000-2802 

2.28 

8.83 

8.98 

<9.04 



3.58 

5.06 

<5.81 

>-23.6 

<10.2 

<10.36 

<9.65 

W2021-2611 

2.44 

9.05 

9.59 

10.23 



5.94 

20.61 

89.96 

>-25.4 

<10.9 

11.24 

10.56 

W2040-3904 

a 

9.00 

9.23 

9.57 



5.30 

8.99 

19.81 

>-24.8 

<10.7 

11.09 

10.41 

W2059-3541 

2.38 

9.23 

9.45 

<9.50 

9.07 

9.10 

8.99 

14.93 

<16.76 

>-25.2 

<10.9 

<10.35 

<9.62 


Redshift assumed to be 2 if no spectroscopic redshift available. 


Mpc scale), double lobes (~7%). Another small percent¬ 
age shows evidence of small scale double lobes on scales 
~2~10 kpc. 

6 .1. The AGN-Heated Source 

In this work we have derived a range in the plausible 
AGN luminosities depending on whether we assume that 
any of the far-infrared luminosity component is AGN- 
heated: 

Abol.AGN = -bAGN.MIR + IaGN^BB 

where Lboi.AGN is the bolometric luminosity of the AGN, 
Lagn.mir is the mid-infrared AGN luminosity from the 
fit to Equation (2), /agn is the fraction of the far-infrared 
emission contributed by the AGN and Lbb is the far- 


infrared luminosity, /agn is assumed to be 0 for the 
Lagn models and to be 1 for the scenario in which we as¬ 
sume that LTotai-Best is Completely AGN-powered. The 
warm AGN component luminosity is very insensitive to 
the dust temperature of the BB component, due to the 
strong constraints placed on the AGN model shape by 
the four WISE data points. The range of values for the 
ratio of the maximum plausible AGN luminosity to the 
minimum estimate is 1.1 < LTotai-Best/LAGN-Min < 5.5, 
with a median value of 2.0. 

We have assumed a covering factor of warm 300K) 
dust that emits in the mid-infrared, 

Tbol.AGN = TAGN.Mm/flwD 

of unity and spherical symmetry, consistent with late- 
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Table 6 

Star Formation Rates and Gas Depletion Timescales 


WISE Name 

Redshift 

30K 

log SFR 
(Le/yr) 

50K 90K 

120K 

log (Gas Depletion Time) 
(years) 

30K 50K 90K 

W0304-3108 

1.54 

1.73 

2.70 

3.16 


9.10 

7.81 

7.04 

W0306-3353 

0.78 

1.22 

2.31 

2.52 


9.29 

7.93 

7.44 

W0354-3308 

1.37 

<1.38 

<2.38 

<2.83 


<9.14 

<7.83 

<7.08 

W0404-2436 

1.26 

1.70 

2.72 

3.17 


9.16 

7.84 

7.09 

W0409-1837 

0.67 


<2.23 

<2.64 



<7.94 

<7.26 

W0417-2816 

0.94 

1.72 

2.79 

3.22 


8.91 

7.56 

6.84 

W0439-3159 

2.82 

2.28 

3.08 

3.57 

4.41 

8.87 

7.66 

6.81 

W0519-0813 

2.05 

<1.54 

<2.44 

<2.92 

<3.82 

<8.99 

<7.74 

<6.93 

W0525-3614 

1.69 

<1.53 

<2.49 

<2.95 

<3.99 

<9.06 

<7.77 

<6.99 

W0526-3225 

1.98 

2.64 

3.55 

4.03 


9.00 

7.74 

6.93 

W0536-2703 

1.79 



3.18 




6.98 

W0549-3739 

1.71 


2.58 

3.05 


9.06 

7.78 

6.99 

W0612-0622 

0.47 

1.13 

2.29 

2.00 


9.40 

7.99 

7.57 

W0613-3407 

2.18 

<1.67 

<2.55 

<3.03 

<3.92 

<8.97 

<7.73 

<6.91 

W0614-0936 

a 

<1.64 

<2.55 

<3.02 

<3.93 

<9.00 

<7.74 

<6.94 

W0630-2120 

1.44 

1.96 

2.95 

3.41 


9.12 

7.82 

7.05 

W0642-2728 

1.34 

1.58 

2.59 



9.14 

7.82 


W0652-2006 

0.60 

1.35 

2.48 



9.35 

7.96 


W0702-2808 

0.94 

<1.34 

<2.41 



<9.24 

<7.89 


W0714-3635 

0.88 

1.42 

2.50 



9.27 

7.91 


W0719-3349 

1.63 


2.99 

3.46 



7.11 

6.28 

W0811-2225 

1.11 



<2.89 




<7.13 

W0823-0624 

1.75 

<1.60 

<2.54 

<3.00 


<9.04 

<7.77 

<6.99 

W1308-3447 

1.65 

1.46 

2.42 

2.89 


9.07 

7.78 

7.00 

W1343-1136 

2.49 


2.67 

3.16 

4.02 


7.69 

6.84 

W1400-2919 

1.67 

<1.02 

<1.97 


<3.40 

<9.07 

<7.79 


W1412-2020 

1.82 

1.76 

2.69 

3.16 


9.04 

7.77 

6.98 

W1434-0235 

1.92 

<1.32 

<2.24 

<2.71 

<3.63 

<9.02 

<7.75 

<6.96 

W1439-3725 

1.19 

<1.03 

<2.05 

<2.50 


<9.18 

<7.86 

<7.11 

W1500-0649 

1.50 


3.08 

3.57 



7.79 

6.99 

W1510-2203 

0.95 

<0.98 

<2.05 

<2.47 


<9.25 

<7.90 

<7.19 

W1513-2210 

2.20 

2.10 

2.98 

3.46 


8.97 

7.72 

6.90 

W1514-3411 

1.09 

<1.05 

<2.10 

<2.53 


<9.20 

<7.86 

<7.14 

W1521+0017 

2.63 



2.87 




8.33 

W1541-1144 

1.58 

1.37 

2.34 

2.81 


9.09 

7.80 

7.02 

W1634-1721 

2.08 

<1.32 

<2.22 

<2.70 

<3.60 

<8.99 

<7.73 

<6.92 

W1641-0548 

1.84 


2.65 

3.12 



7.76 

6.97 

W1653-0102 

2.02 

<1.28 

<2.18 

<2.65 

<3.55 

<9.00 

<7.75 

<6.95 

W1657-1740 

a 

<1.27 

<2.18 

<2.66 

<3.56 

<9.01 

<7.75 

<6.94 

W1702-0811 

2.85 

<1.52 

<2.31 


<3.64 

<8.86 

<7.66 


W1703-0517 

1.80 

1.36 

2.29 

2.72 

3.69 

8.56 

7.23 

6.44 

W1707-0939 

a 

<1.39 

<2.30 

<2.78 

<3.67 

<9.00 

<7.74 

<6.93 

W1936-3354 

2.24 

1.69 

2.56 


3.93 

9.28 

8.08 


W1951-0420 

1.58 

<1.32 

<2.29 


>5.41 

<9.08 

<7.79 


W1958-0746 

1.80 

<1.32 

<2.25 

<2.72 

<3.65 

<9.04 

<7.77 

<6.98 

W2000-2802 

2.28 

<0.96 

<1.83 

<2.31 

<3.20 

<9.40 

<8.16 

<7.34 

W2021-2611 

2.44 

2.09 

2.92 


4.28 

9.15 

7.97 

0.00 

W2040-3904 

a 

2.09 

3.00 

3.47 


9.00 

6.94 


W2059-3541 

2.38 



<2.70 

<3.11 



<6.92 


2 = 2 assumed for sources without a spectroscopic redshift. 


stage mergers with heavily-obscured nuclei. This may 
not be a valid assumption in which case we may have 
under-estimated the AGN luminosities and BH masses. 
In particular narrow emission lines photoionized by the 
AGN are seen in many cases, therefore these luminosity 
estimators may be excluding flux emitted by the AGN 
that does not intercept the dusty structure. A further 
complication is that the covering factor of the cold 
50K) dust is likely to be different from that of the hot 
dust, and its optical depth to the mid-infrared emission 
is unknown. 

An added complication is the possibility that a non- 
spherical source emits non-isotropically because it is op¬ 
tically thick, as is often the case for torus models (eg. 
Efstathiou et al. 2013). In that case a further correc¬ 
tion for anisotropy is required, and this depends on the 


particular torus model and also on the inclination of the 
torus. Generally speaking, for edge-on viewing the lumi¬ 
nosity will be underestimated whereas for face-on view¬ 
ing it may be overestimated. Face-on viewing is ruled 
out by the red MIR-optical SEDs and therefore we can 
conclude that Lagn.mir provides a firm lower limit to 
the bolometric luminosity of the AGN. 

Support for the high covering factor interpretation 
for most of our WISE-NVSS-ALMA sample comes from 
the overall relative numbers of obscured and unobscured 
AGN amongst the highest luminosity radiative mode 
AGN at high redshift (z^ 1-3), which are roughly equa l 
(ILacy et al. 2015 Assef et al. 2015 Tsai et al. 2015). 
Taken as face value this would imply an average cover¬ 
ing fraction of ~50%, if these populations differ only by 
orientation (ignoring the differing selection functions for 
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type 1 and type 2 AGN, Elitzur 2012). 

6.1.1. A Torus? 

Although we have commented above on the likelihood 
that the covering factor may be lower than 1.0 for many 
sources, it might be questionable to conclude that a clas¬ 
sical smooth torus could be responsible for the very high 
luminosities found for some sources, because they would 
require very large tori (several hundred pc to over a kpc 
in diameter), especially if the Axotai-Best values are in¬ 
terpreted as fully AGN-powered. Such large thin struc¬ 
tures would be unstable. A clumpy torus would be more 
plausible, as it could achieve a wider range of dust tem¬ 
peratures than a smooth torus of the same diameter. 

A torus-like structure would imply that a large frac¬ 
tion of the total AGN luminosity escapes dust absorption 
and would be easily visible at optical wavelengths. Since 
the line-of-sight optical continuum emission is very faint 
compared to type 1 AGN (Figure the tori would all 
have to be inclined closely to our line of sight. Our selec¬ 
tion function of course has favored the selection of highly 
obscured systems, and we do indeed find a fraction of our 
sample to have (s mall) double radio lobes in the plane 
of the sky (Section 5.4), as might be expected for a ra¬ 
dio quasar or radio galaxy interpreted as viewed through 
an edge-on disk in the standard unification picture. The 
majority of our sample does not display extended radio 
lobes however. 

6.1.2. Extended NLRs, Polar Winds and Lyman-a Blobs 

High resolution observations of nearby low luminos¬ 
ity AGN show that the MIR emission can lie in the po¬ 
lar direction and may be associated with NLR clouds or 
with thermal winds from the AGN (Zhang et al. 2013 
and references therein). It is also known that a signif¬ 
icant fraction of FRII radio AGN have extended e mis¬ 
sion line regions (EELRs) up to tens of kpc in size (Fu 
fc Stockt ^ 20071), which are outflowing (|Shih & btock- 
ton|2014[ ). I'he Gompact-Steep-Spectrum (GSS) sources, 
which are younger versions of the FR IIs, have smaller 
EELRs which are bet ter aligned with t he radio struc¬ 
tures than in FR IIs (|Axon et al. 2000). Therefore it 
is possible that some of the iVllR emission in our sam¬ 
ple is associated with NLR or EELR clouds which are 
heated by the centr al AGN, or shock-heat ed by the ra¬ 
dio je t inter actions (Mullaney et al. 2013). Efstathiou 
et al. (2013) have found that a dusty MLR component 
is needed to fit the MIR-FIR SED of the hyperluminous 
g alaxy IRAS10214-H47 24 at z = 2.285. 


Bridge et al. (2013) have discovered that a subset of 
the Hot UUGs which have extended faint optical emis¬ 
sion also possess Lyman Alpha Blobs (LABs; defined to 
have Lya emission extended on scales > 30 kpc). The 
WISE LABs are non-symmetric and there is evidence 
for large outflow velocities. The presence of extended 
ionized gas suggests significant shock heating or that a 
significant fraction of the nuclear ionizing radiation must 
be able to find its way out of the galaxy, as in the ra¬ 
dio galaxy EELRs, even in these highly obscured MIR- 
dominated systems. Few of our sources currently have 
spectra that cover Lya. W0613-3407 has Lya extended 
spatially on a scale of 3", so qualifies as an LAB by 
the size definition. W1343-1136 appears to also have ex¬ 
tended emission (about 2.5"; 25 kpc) but to just miss 


the usual LAB definition of >30 kpc. Although it re¬ 
mains to be seen whether a large fraction of our WISE- 
NVSS-ALMA sources possess LABs, we have some evi¬ 
dence that some of them possess broad forbidden emis¬ 
sion lines th at might indicate substantial outflows. |Kim| 


et al. 
broac 


(2013) find that the [O III] lines are exceptionally 
or SIX of our quasars, with full width at half max- 
•^1300 to 2100 km s“^, significantly larger than 


imum 

that of typical distant quasars. 


6 .2. Star Formation 

The contribution to the 1-lOOO^m luminosity from star 
formation has larger overall uncertainty than the AGN 
contribution because star formation can produce a larger 
range in observed dust temperature, varying by ^2 dex 
between the viable 30 K to 120 K BB models (Table 
columns 4-6). The maximum star formation luminosity, 
LBB-Max, is comparable to the the AGN luminosity esti¬ 
mators, while the minimum starburst luminosity, coming 
from the 30 K model in all but one case, can be well over 
an order of magnitude smaller than the AGN luminos¬ 
ity. A minimum SFR of 0 is also possible if Lxotai-Best 
is interpreted as completely AGN-powered. 

If we adopt 50K dust for the BB component, the SFRs 
lie between ~200 and 3500 Mq yr“^. The lower val¬ 
ues are consistent with main sequence galaxies at these 
redshifts while the high er values may require a starburst 
( Delvecchio et alT]|20i5 ). 


Uur radiative transfer modeling of our best observed 
source favored the presence of a young (10-15 Myr), pre¬ 
sumably compact, starburst contributing about 50% of 
the 1-1000/im luminosity. Gompact starbursts have been 
investigated in several local galaxies. The nuclear struc¬ 
tures around the twin nuclei of Arp 220 are thought to 
be powered by compact starbursts with relatively high 


callv thick out to >>100 um (Wilson et al. 

2014 Scov- 

ille et al. 

2015 

Barcos-iVlunoz et ai. 


'I'sai et al. 


]2015) have considered in more detail the possibility that 


a compact starburst could contribute a significant lumi¬ 
nosity to their most luminous Hot DOGs, those with 
Lboi > lO^'^L©. Using He 2-10 as a local analo g and the 
STARBURST99 code (Leitherer & Ghen 2009), they find 
SFRs > 5xlO^“^M0 yr“^ for the conservative case of a 
t op heavy IMF. 


Tsai et al. (2015) conclude that there is insufficient GO 
in these systems to support very large amounts of star 
formation. In our case we find mean ISM masses of be¬ 
tween ~ 2-9 xIO^^Mq, depending on the assumed dust 
temperature. For the range of estimated SFRs across 
our models the gas would be depleted in ~ 2 Myr - 2 
Gyr. The lower depletion times (corresponding to the 
higher SFRs) are not insupportable, in the scenario of a 
late-stage, violent, gas-rich merger. Therefore it is quite 
possible that a vigorous starburst is present in some of 
our systems, and that a compact young starburst may 
contribute to the warmest dust emission. 

6.3. High Accretion Rates Relative to Star Formation 

In Figure we plot LAGN-Min vs. Lbb-sok: no corre¬ 
lation is apparent. We compare our sample t o the radio¬ 
quiet 6 'pztzer-IRAG-selected quasar sample of Lacy et al. 
(2011) (Lacyll), for which fits have been done using the 
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Figure 9. Comparison of the Spftzer-selected sample of Lacyll 
with our sample in the Lagn vs . Lbb plane. For this figure we 
show the BB luminosity from our 50K model to match the method 
used by Lacyll. 

same formalism as here. The Lacyll sample has a wider 
range in mid-IR color selection than our sample, includ¬ 
ing Type 1 quasars, reddened Type 1 quasars and Type 2 
quasars, and is therefore representative of the IR-bright 
AGN population as a whole. There is some evidence that 
the heavily obscured quasars in our sample have similar 
far-infrared luminosities to the Lacyll sample, ie. simi¬ 
lar star formation rates, but systematically higher AGN 
luminosity. This is consistent with them being in a sys¬ 
tematically higher accretion rate phase relative to star 
formation. 


6.4. ISM Masses 

The masses of the BH, the ISM and the stellar com¬ 
ponent are tabulated in Table The ISM mass is di¬ 
rectly proportional to the 870 /im luminosity, but de¬ 
pends much less strongly on dust temperature than does 
Lbb and SFR. 

We found in Section [ST] that the ISM masses are com- 


in the GOSMOS field ( 

Scoville et al. 

2014). We can 

also c 
molec 
(2013 

ompare our sample to the compila 
ular masses for z>l systems of[ 
), their Figure 9. Assuming a GO 

tion ot all known 
(Jarilli & Walter] 
to H 2 conversion 


1998 Bolatto et al. 2013) of aco~ 0.8 iVl© / (K km 


s ^ pc^), and a ssuming 50% of the g aseous ISM is in 


molecular form (Garilli & Walter!2013), our median de¬ 


tected ISM mass implies a median molecular gas mass 
of 2.95xlO^°M0 and a median GO line luminosity of 3.7 
xl0^° (K km pc^). This lies at about the 80*^ per- 
centile of t he hig h redshift QSOs in Figure 9 of |Garilli| 
& Walter (|2013), and about the 50*^ percentile of the 
SiVIGs. It is lower than four of the RGs in this figure and 
comparable to the two others. 

In summary, if we assume no contribution to the mea¬ 
sured ALMA fluxes from non-thermal synchrotron emis¬ 


sion, we find that the implied ISM masses are quite high, 
comparable to those of typical “IR Bright” star forming 
systems in the GOSMOS field at these redshifts, to the 
most gas ric h galaxies in the local GO survey of |Leroy| 
et al. (2009|), and consistent with the large CO masses 
ot z>l OSOs, SMGs and RGs. This gas and dust could 
exist in nuclear, AGN-heated structures, or it could be 
powered by star formation somewhere within the host 
system. 


6.5. BH Masses 


The BH mass estimates derived by Kim et al. (2013) 


for 3 of the 6 quasars with available [O III]A5007A line 
luminosities are significantly larger that the values de¬ 
rived from our MIR data. Our masses may underesti¬ 
mated due to extreme extinction in the MIR, (although 
we might expect to recover such dust-absorbed energy in 
the FIR-submm), and our assumption of ~ 100% cover¬ 
ing factor may be incorrec t. Another possi bility is that 
the [O III] line strengths of Kim et al. (2013) are boosted 
by shocks and outflows. We will address the relationship 
between Lmir and L[oiii] for these systems (cf. Mul- 
laney et al. 2013) in a forthcoming paper presenting the 
spectroscopy. 

6.6. Nature of the Radio Sources 

We will address the morphology of the radio sources 
in our paper presenting the VLA results, where we will 
show that the majority of the sample are compact on ^1- 
3 kpc scales. The rest frame 3 GHz radio powers lie in the 
range log (Fs.oghz /WHz~^) = 24.74 - 27.33. The radio 
power of both RQ and RL systems evolves with redshift 
(Best et al. 2014) and some radio sources at z =1-3 with 
radio powers in this range are found to radio qu iet based 
on the 24 u rn q value, ^24 = log(/ 24 Mm// 20 cm) ( [Simpson 
et al. 2012). Our sources have values of <722 that are 
signiticantiy too low for them to be considered radio¬ 
quiet by this crit erion. 

In Figure \T0\-right we compare LAGN-Min with 
Tkadio-sGHz, and also include the hig h-redshift (1 < 
z <5.2) r adio galaxy (HzRG) sample of|De Breuck et 
^ (|2010| , and the 3G radio galaxy sample ot|Gieary et 
aE (2007^ at 0.4< z <1.2. For the HzRG sample we also 
include the estima ted core fiux, using the 20 GHz core 
fractions given by De Breuck et al. (2010) and assum¬ 
ing the same fraction is appropriate at 3 GHz. There is 
a large disparity between our sample, the Lacyll sam¬ 
ple and the two high redshift samples. The HzRGs and 
3G RGs have 2-3 orders of magnitude more radio power 
than our sample, for a given mid-IR AGN luminosity, 
which is expected given that the power of the HzRGs is 
dominated by the extended lobes. The core radio pow¬ 
ers for the classical HzRGs have a similar range as our 
sample but our sources are significantly more luminous 
in the infrared, consistent with a higher accretion rate. 
The radio-blind MIR-selected sample of Lacyll has on 
average 2-3 orders of magnitude less radio power than 
our sample, as expected for an RQ-dominated sample. 

All four samples in this figure display an apparent 
correlation between Lagn and FRadio-sGHz- Given the 
rarity of these sources and the large redshift range in¬ 
volved, the correlation for our sample may be an arti¬ 
fact of Malmquist bias resulting from the 22 gm and 
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Figure 10. Left: Selection criteria for the ALMA quasar sample at 22 and 20 cm. An apparent weak correlation between 22 ^m and 
20 cm flux density may be due to the selection boundaries (dashed lines). Also shown is the 870 ^m flux density vs. 20 cm flux density, 
showing no trend. Right: AGN luminosity vs. r adio power for our W ISE sample and for the HzRGs of|De Breuck et al.|l|2010ll (both total 
emission and core emission) the 3G quasars from lCleary et ar] ( |2007| l and the MIR-selected AGN of Lacyli. 


the 1.4 GHz flux density thresholds. We show in Fig¬ 
ure ^left that there exists an apparent correlation 
between 22 and 20 cm flux density that may be 
caused by the limited dynamic range in the flux ratio 
selection thresholds for the WISE-NVSS-ALMA sam¬ 
ple: -1 < ^ 05 (/ 22 Mm// 20 cm) < 0 (dashed lines). There¬ 
fore the apparent cor rela tion seen between Lagn and 
-fRadio- 3 GHz in Figuremay not be real, at least 
for our sample. 


7. CONCLUSIONS 

We have selected a sample of extremely red, luminous, 
radio-powerful sources in the 0.5 < z < 3 redshift range 
using WISE MIR colors and the NVSS & FIRST 20 cm 
radio surveys. We present ALMA 870 /rm photometry 
for 49 southern sources from the total sample of 156 red 
sources, and redshifts for 45 of them from a combination 
of optical and near-infrared spectroscopy. JCMT imag¬ 
ing at 850 /xm has been pres ented for 30 of the northern 
sample by Jones et al. (2015). Combined, 30 sources have 
a detection at 850 or 870 gm, 25 of them with a known 
redshift. We also have i?-band imaging for 27 sources, 
Herschel photometry for fifteen sources, including two 
of the ALMA sample, and CSO 350 /xm data four addi¬ 
tional sources, including two from the ALMA sample. 

Having compared the SEDs of the red WISE selected 
sources with other samples and template SEDs, we con¬ 
clude that the rest frame MIR-submm SEDs of our 
WISE-NVSS sources are dominated by AGN emission 
in the MIR. They have extremely red optical-MIR col¬ 
ors and high bolometric luminosities in the ULIRG and 
HyLIRG regime. In some cases approaching or exceeding 
101 ^ 0 . They are redder in the NIR-MIR than previous 
samples selected from Spitzer surveys, including almost 
all of the Spitzer DOGs. They also display system¬ 
atically warmer overall MIR-submm SEDs, and proba¬ 
bly have higher levels of accretion, relative to star for¬ 


mation, than the DOGs. BH mass estimates for our 
sample are 7.7 < log(MBH/Al 0 ) < 10.2. We conclude 
that these sources are best labeled as obscured radio¬ 
powerful QSOs. The rest-frame 3 GHz radio powers are 
24.7 < log(P 3 .oGHz/IFi 7 z“^) < 27.3, and all sources are 
radio-intermediate or radio-loud. The ability of WISE 
to find this rare and unique sample is due to the large 
volume accessible to WISE. 

Our best constrained source has radiative transfer so¬ 
lutions with ^ equal contributions for an obscured AGN 
and a young (10-15 Myr) compact starburst. Simpler 
two-component fits to the whole sample And that the 
SFRs of the sample could be in the hundreds to thou¬ 
sands of solar masses per year range, but it also possible 
to fit a significant fraction of the ALMA data for some 
sources with a centrally-heated dusty structure, in which 
case most of the entire bolometric luminosity could be 
attributed to the obscured AGN. In that scenario it is 
likely that the emission is dominated by a small high 
covering factor cocoon and/or an extended NLR. 

Our sample is simila r in MIR se l ection met hod to 


Eisenhardt et al. (2012), Wu et al. (2012) and Bridge 


et al. 12013), who did not use radio flux density as a 


selection criterion. Their Hot DOG samples exhibit sim¬ 
ilar SED shapes to our radio-selected sample. The MIR 
SEDs of the EWB12 samples may be steeper on average 
than those of our radio-selected sample, however this may 
be due to different selection effects between the MIR- 
optical color selection criteria for the two samples, re¬ 
sulting in an average higher redshift for the radio-blind 
sample. 


This paper makes use of the following ALMA data: 
ADS/JAO.ALMA#2011.0.00397.S. ALMA is a partner¬ 
ship of ESO (representing its member states), NSF 
(USA) and NINS (Japan), together with NRG (Canada), 
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NSC and ASIAA (Taiwan), and KASI (Republic of 
Korea), in cooperation with the Republic of Chile. 
The Joint ALMA Observatory is operated by ESO, 
AUI/NRAO and NAOJ. The National Radio Astron¬ 
omy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by 
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use of data products from the Wide-held Infrared Sur¬ 
vey Explorer, which is a joint project of the Univer¬ 
sity of California, Los Angeles, and the Jet Propulsion 
Laboratory/California Institute of Technology, funded 
by the National Aeronautics and Space Administra¬ 
tion. This work is based on observations made with 
the Caltech Submillimeter Observatory, which is oper¬ 
ated by the California Institute of Technology under 
funding from the National Science Foundation, contract 
AST 90-15755. This paper uses data from SDSS (DR 
8 ). Funding for SDSS-III has been provided by the Al¬ 
fred P. Sloan Foundation, the Participating Institutions, 
the National Science Foundation, and the U.S. Depart¬ 
ment of Energy Office of Science. The SDSS-III web 
site is http://www.sdss3.org/, RJA was supported by 
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prove the paper. 
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Figure 8. continued. SED plots. 
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Figure 8. continued. SED plots. 
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Figure 8. continued. SED plots. 
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Figure 8. continued. SED plots. 
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Figure 8. continued. SED plots. 
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Figure 8. continued. SED plots. 












